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1 Introduction 

The inverse spectral problem on a Riemannian manifold (M, g) , possibly with 
boundary, is to determine as much as possible of the geometry of (M, g) from 
the spectrum of its Laplacian Ag (with some given boundary conditions) . The 
special inverse problem of Kac is to determine a Euclidean domain Q C M" up 
to isometry from the spectrum Specsi^) of its Laplacian with Dirichlet, 
Neumann or more general boundary conditions B. The physical motivation 
is to identify physical objects from the light or sound they emit, which may 
be all that is observable of remote objects such as stars or atoms. 

The inverse spectral problem is just one among many kinds of inverse 
problems whose goal is to determine a metric, domain or scatterer from phys- 
ically relevant invariants. A comparison with other inverse problems shows 
just how small a set of invariants the spectrum is. For instance, the boundary 
inverse problem asks to determine the metric on a fixed bounded domain 
C M of a Riemannian manifold (M, g) from the spectrum of the Dirichlet 
Laplacian on L'^{Q), and from the Cauchy data ^^\dn of its eigenfunctions, 
or equivalently from its Dirichlet-to-Neumann operator [BK, KKL, LU]. The 
inverse scattering problem seeks to determine an obstacle from its scattering 
amplitude [Ma, Ma2], or from its scattering length spectrum [St, St2], a set 
of lengths parametrized by 5"""^ x S'^~^. By comparison, the inverse spectral 
problem has to make do with just a discrete 'unformatted' set of eigenvalues 
(or resonance poles in the open case). 

*Research partially supported by NSF grants # DMS-0071358, DMS-0302518. 



2 



It is probably a consequence of this relative poverty of invariants that 
most of the results in inverse spectral theory over the last two decades (es- 
pecially since the appearance of Sunada's aricle [Su]) are 'negative results', 
showing that one cannot determine metrics or domains by their spectra. The 
collection of non-isometric isospectral pairs of Riemannian manifolds would 
require a lengthy survey of its own and for that we refer to the lectures of C. 
Gordon (cf. [Gor, Gor2]). By comparison, the number of 'positive results' 
showing that one can indeed recover a domain or metric is rather small. 

This survey is devoted to the positive results. The emphasis is on relations 
between the spectrum of the Laplacian and the dynamics of the geodesic flow 
G* : S*M — > S*M. Most of the new material concerns wave trace invariants 
and their applications to solving concrete inverse spectral problems. The 
wave group is the quantization of the geodesic flow, and so wave trace meth- 
ods often 'reduce' inverse spectral problems to inverse dynamical problems. 
Because of their relevance, we have attempted to describe inverse dynamical 
problems and results. 

1.1 Some basic inverse spectral problems 

Let us introduce some basic terminology. The spectrum of a compact Rie- 
mannian manifold defines a map 



from some class of metrics M on a manifold M to the spectrum of its Lapla- 
cian, 



with boundary conditions B : C°°{M) C^{dM) if dM ^ 0. Here, A 
denotes the positive Laplacian 



of a Riemannian manifold {M,g), where gij = g{-^., [5'*-'] is the inverse 
matrix to [gij\ and g — deffg^y]. We will only consider Dirichlet Bu — u\qm 



Spec : M ^ M^, {g, B) ^ ^pec(A,,s) = {Ao < A? < < • • • } 
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and Neumann Bu — di,u\QM- Eigenvalues are repeated according to their 
multiplicities. 

Two metrics or domains are called isospectral if they have the same spec- 
trum. The main problem in inverse spectral theory is to describe the possible 
spectra A C of Laplacians and, for each possible spectrum, to describe 
the metrics or domains in the spectral class 

Spec-\A). (1.1) 

Somewhat simpler is to describe the possible smooth curves in the isospectral 
class, since it apriori eliminates irregular subsets. An isospectral deformation 
of a Riemannian manifold (possibly with boundary) is one-parameter family 
of metrics satisfying Spec{M, gt) = Spec{M, qq) for each t. Similarly, an 
isospectral deformation of a domain with a fixed background metric qq and 
boundary conditions S is a family with SpecsiP't) = SpecBiyt)- One 
could also pose the inverse spectral problems for boundary conditions (while 
holding the other data fixed) as in [GM2, PT]. 

The inverse spectral and isospectral deformation problems are difficult 
because the map Spec is highly nonlinear. The linearization of the problem 
is to find infinitesimal isospectral deformations, i.e. deformations for which 
the eigenvalue variations vanish to first order. By first order perturbation 
theory, the variations of the eigenvalues under a variation of the metric are 
given by 

>^3 = Jt^Mt=o = (A0,-,0,), (1.2) 

where A is the variation of the Laplacian and where (pj = 4>j{0) is an or- 
thonormal basis of eigenfunctions which varies smoothly in t (such a basis 
exists by the Kato-Rellich theory). As will be recalled below, lengths of 
closed geodesies 7 are spectral invariants (at least, when there is at most one 
closed geodesic of each length), so infinitesimal iso-length-spectral deforma- 
tions are those where 

= 0, V7. (1.3) 

The deformation of the metric is a symmetric 2-tensor g, so the linearized 
problem is to determine the space oi g & S^T*M (modulo tensors arising 
from diffeomorphisms (t>l{g)) for which 

fgds = OW^^ = 0, (Vj). (1.4) 

J 'J 
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The operator A = Op{g) is the difTerential operator with symbol g. The 
hnearized problem is still very difficult because it requires a study of the 
asymptotic behavior of the expressions (1.4) as the lengths or eigenvalues 
tend to infinity. This is tantamount to the study of the equidistribution 
theory of closed geodesies and eigenf unctions. 

The basic distinctions in inverse spectral theory are the following. We 
say that 

• a metric or domain is spectrally determined (within M) if it is the unique 
element of M with its spectrum; 

• it is locally spectrally determined if there exists a neighborhood of the 
metric or domain in M on which it is spectrally determined; 

• a metric or domain is spectrally rigid in M if it does not admit an 
isospectral deformation within the class; 

• the inverse spectral problem is solvable in M if Spec\y[ is 1 — 1, i.e. if 
any other metric or domain in M with the same spectrum is isometric 
to it. If not, one has found a counterexample. 

There are analogous problems for Laplacians on non-compact Riemannian 
manifolds, which often have continuous spectra as well as discretely occurring 
eigenvalues. In place of eigenvalues, one considers the resonances i?es(A) of 
A, i.e. the poles of the analytic continuation of its resolvent 

R{z) = (A + ^2)-i_ 

Depending on whether the dimension is odd or even, i?es(A) is a discrete 
subset of C or of the logarithmic plane. The inverse spectral problems above 
have natural analogues for resonance poles. We refer to Zworski's expository 
articles [Zw2, Zw3] for background. 

To illustrate the current state of knowledge, we note that even the sim- 
plest special metrics are not known to be spectrally determined at the present 
time (at least to the author's knowledge). It is not known: 

• if the standard metric g^ on 5"" is determined by its spectrum (in di- 
mensions > 7), i.e. if {M,g) (or even {S",g)) is isospectral to {S'",go) 
then it is isometric to it. This has been proved in dimensions < 6 [T]. 
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• if ellipses in the plane are determined by their Dirichlet spectra, or even 
if they are spectrally rigid, i.e. if there exist isospectral deformations 
of ellipses (with Dirichlet boundary conditions) . 

• if hyperbolic manifolds are determined by their spectra in dimensions 
> 3. I.e. if (Mo,5'o) is hyperbolic and {M,g) is isospectral to it, then 
is {M,g) hyperbolic? This is of course true in dimension 2. Is {M,g) 
isometric to {Mo,go)7 In dimension 2, this is known to be false for 
some hyperbolic surfaces. 

• if flat metrics are determined by their spectra in the sense that if (M, go) 
is flat and (M, g) is isospectral to it, then (M, g) is flat (it is known 
that this is true in dimensions < 6 or in all dimensions if additionally 
g is assumed to lie in a sufficiently small neighborhood of go [Ku3]); 
it is also classical that there are non-isometric flat tori with the same 
spectra. 

These special cases are tests of the strength of the known methods. An- 
other test is given by the two-dimensional inverse spectral problem. One- 
dimensional problems are comparatively well understood because the eigen- 
value problems are ordinary differential equations and the underlying dy- 
namics consists of just one orbit (an interval)! Two-dimensional problems 
are already rich in spectral and dynamical complexities, as illustrated by the 
classical dynamics of twist maps or geodesic flows on Riemannian surfaces. 
In general, the inverse spectral problem grows rapidly in difficulty with the 
dimension, and is already quite open for analytic surfaces and domains in 
two dimensions. This motivates our concentration on two-dimensional prob- 
lems for much of the survey. The following simple-sounding problems are 
still apparently beyond the reach of known methods: 

• Are convex analytic domains determined by their spectra among other 
such domains? Are they spectrally rigid? 

• Are convex analytic surfaces of revohition determined by their spectra 
among all metrics on 5*^? Are they spectrally rigid? 

These problems are in some ways analogous to each other in that the 
unknown is a function of one variable (the boundary or the profile curve), 
and that is limit of what wave trace invariants at one orbit can hope to 
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recover. Surfaces of revolution are simpler than plane domains since the 
geodesic flow is integrable, while billiards on plane domains could have any 
dynamical type. On the other hand, in the domain problem the class M only 
consists of convex plane domains, whereas in the second problem we allow 
any other metric. If we similarly restricted the class of surfaces of the first 
problem entirely to analytic convex surfaces of revolution, then the answer 
is known to be 'yes' [Z2]. 

The inverse problem for surfaces of revolution might sound reasonably 
simple since the geodesic flow is completely integrable, and the feeling arises 
that one should be able to detect this property from the spectrum. This is 
one of many problems which relate spectral invariants to dynamics of the 
geodesic flow. 

The relations between Laplace spectrum and dynamics have been at the 
center of at least the positive results in inverse spectral theory in the last 
thirty years, by comparison with the emphasis on heat invariants in the 
earher period. Before going into the technical relations between spectral (in 
particular, wave) invariants and dynamics, it might be helpful to give some 
heuristic principles which suggest the relevance of dynamical inverse spectral 
problems to Laplace inverse spectral problems. The first is the relation of 
classical to quantum mechanics. Two Laplacians are isospectral if 

A,, = C/A,,C/*, (L5) 

where U : L'^{Mi,gi) — > L'^{M2,g2) is a unitary operator. In the Dirac dic- 
tionary of analogies, the classical analogue of this similarity is the symplectic 
conjugacy 

of the corresponding geodesic flows (cf. (3.8)). Here, x • T*Mi\0 T*M2\0 
is a homogeneous symplectic diffeomorphism.) This analogy should not be 
taken too literally, but it is useful in suggesting conjectures. In modern 
language, the analogue would hold if U were a unitary Fourier integral op- 
erator quantizing x- K would clearly be difficult to prove, even in special 
cases, that isospectral Laplacians are conjugate by unitary Fourier integral 
operators, though it was observed independently by Uribe and the author 
(see [Z6]) and by P. Berard ([Be, Be2]) that the Sunada counterexamples 
[Su] have this property (the resulting Fourier integral operators were termed 
'translplantations' by Berard). It was also observed in [Z6] that such Fourier 
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integral intertwining operators need not be quantizations of symplectic dif- 
feomorphisms, but could be (and indeed they are, in the Sunada examples) 
quantizations of multi- valued symplectic correspondences. 

But the analogy is suggestive and is fruitful on a local (or more accurately, 
formal local) level. One of the main results in inverse spectral theory in recent 
years is the theorem due to V. Guillemin [G, G3] (see also [Z3, Z4]) that 
isospectrality (with a simple length spectrum assumption) implies the formal 
local symplectic equivalence of the geodesic flows around corresponding pairs 
of closed geodesies, i.e. it implies the equality of their Birkhoff normal forms. 
This implies local symplectic equivalence around hyperbolic orbits, although 
not around elliptic orbits (see Problem 3.4). It is perhaps the closest that 
the above heuristic principle comes to being valid in a general setting. 

A further heuristic principal is that much of the dynamics of a flow is 
encoded in the structure of the flow near closed orbits. To the extent that 
this is true, local equivalence would be a powerful fact and one should be able 
to obtain strong information about the metric g from studying the wave trace 
expansion around closed geodesies. This raises the possibility that the zeta 
functions of the flows might determine their dynamical type (see Problem 
3.5). 

• Can one determine the dynamical type of the geodesic flow from the 
spectrum of A, i.e. whether the geodesic flow is integrable, ergodic or 
of some other type? The Rucllc zeta function of the geodesic flow is 
generically a Laplace spectral invariant. Can one determine dynamical 
type from its analytic properties. 

• In the above problems, we almost always assume in addition that the 
length spectrum is multiphcity-free: i.e. that the set of closed geodesies 
of a fixed length, or more generally the fixed point sets of the geodesic 
fiow, should contain at most two components interchanged by the time 
reversal involution (x, ^) — * (x, — ^). Are any of the known counterex- 
amples, i.e. non-isometric isospectral pairs, multiplicity free? 

The motivation for the second problem is that wave trace methods cannot 
get off the ground, and in particular the dynamical zeta function need not 
be a spectral invariant, unless the length spectrum is simple. The issue is 
that there could exist complicated cancellations among invariants of closed 
geodesies of the same length. 
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1.2 Strategies for solving the inverse spectral problem 



As mentioned above, we are concerned here primarily with the positive re- 
sults, ones which prove that certain geometric data is determined by the 
spectra. Roughly speaking, the strategy for obtaining positive results has 
long consisted of the the following steps: 

(A) Define a lot of spectral invariants; 

(B) Calculate them in terms of geometric or dynamical invariants; 

(C) Try to determine the metric or domain from the invariants. 

A crucial limitation arises in step (B), which accounts for the relative 
paucity of positive results compared to negative results. It is easy to define 
a complete set of spectral invariants, namely the 'special values' of any one 
of 

' The heat trace, Z(t) = Tre"*^ = J2T=o > 0), 

< The zeta function = ^^rA"^ = ^JIq Aj^" (3fJs > n) 

The wave trace S{t) = Tre'^^ = ^J^q e*^^*, or Sev{t) = Tr cost^A. 

(1.6) 

Of course, C(s) must be meromorphically continued to C and S{t) is a 
distribution rather than a function. But the point we are making is that spe- 
cial values are rarely computable in terms of the geometry and are therefore 
of limited use for positive results on the inverse spectral problem. By com- 
parison, special values can be used to prove negative results by showing that 
the traces of any of the above operators are the same for two non-isometric 
(M,g). A key step in obtaining positive results is to find computable in- 
variants, and to give efficient algorithms for computing them in terms of the 
simplest possible geometric invariants. It should be mentioned that there 
exists spectral invariants such as Ai and logdet A = — C'(0) which are useful 
although not computable in the above sense. To maintain our theme of wave 
invariants and dynamics, we will not discuss such invariants. 

Oversimplifying a bit, the computable invariants arise from the singularity 
(or related) asymptotic expansions of the traces defined above or, in another 
language, from non-commutative residues of functions of the Laplacian (we 
re refer to [G, G2, Z9] for discussion of non-commutative residues). In fact, 
all of the computable invariants known to the author are wave invariants, i.e. 
arising from the singularities of the distribution trace or residues of the wave 
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operator U (t) — e**"^ at times t in the length spectrum of (M, g) (including 
t = 0). 

And it should not be forgotten that the goal of inverse spectral theory is 
step (C). There now exist a number of rather abstract results showing that 
various dynamical or quantum mechanical invariants (e.g. Birkhoff normal 
forms) are spectral invariants [G, GM2, ISZ, Zl, Z3, Z4]. But this only 
trades one inverse problem for another, and there is relatively little work on 
the subsequent inverse problem of determining the domain or metric from 
these invariants. For instance, it is not hard to see that the classical Birkhoff 
normal form of the Poincare map of a bouncing ball orbit does not determine 
all of the Taylor coefficients of the boundary at the endpoints of the orbit, 
even if the domain has one symmetry. Often step (C) is the deepest, requiring 
a separate study of inverse dynamical problems. 

1.3 Contents of the survey 

The focus of this survey is on the use of wave trace formula to derive infor- 
mation about the metric and geodesic flow around closed geodesies and to 
determine metrics or domains from the information. We survey in some detail 
the relation between wave invariants and Birkhoff normal form invariants on 
general Riemannian manifolds, initiated by V. Guillemin, and developed by 
the author and by lantchenko-Sjostrand-Zworski [G, G3, ISZ, SjZ, Z3, Z4]. 
In an appendix to this article [SjZ2], J. Sjostrand and M. Zworski describe in 
more detail how their general results on quantum monodromy apply to the 
Laplacian. 

But as mentioned above, to succeed with step (G) we need to be able to 
determine a metric or domain from such invariants. In the end, the crucial 
problem is to compute wave trace invariants in the simplest and most efficient 
way, to analyze them in detail, and to reconstruct the domain or metric. This 
is most feasible in dimension two, so we review in some detail the articles 
which have succeeded in determining special families of domains or metrics 
from wave invariants, to wit, bounded plane domains [CdV, P, Zl, Z2, Z5, 
Z7, ZIO, ISZ, GM, MM, SI] and surfaces of revolution [Z2]. One of our aims 
is to describe a new method for calculating wave invariants from [Z5] which 
so far has achieved better results than the Birkhoff normal form approach. 

In addition, we provide a fair amount of background that hopefully puts 
the special problems in context. There already exist a number of surveys on 
the inverse spectral problem (e.g. [Ber, Be3, G3, Me, Gor, Gor2]) including 
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our own expository articles [Z9, ZIO], and we have tried to avoid duplication 
of material which already appears elsewhere. However, to make the survey 
more self-contained we quote from a number or prior surveys, including our 
own. We also follow the lecture notes of Melrose [Me] in our discussion of the 
Lifshits (Penrose mushroom) example of two domains with the same wave 
invariants, and also the (much better known) examples of domains with the 
same heat invariants. 

We also omit a number of topics as being too far from our focus on wave 
invariants and dynamics. As mentioned above, we do not discuss counterex- 
amples and negative results (cf. [Gor, Gor2]). We also omit discussion of 
compactness results of isospectral sets, of which there are many since the (un- 
published) work of Melrose and the work of Osgood-Phillips-Sarnak [OPS] on 
isospectral sets of plane domains (for the resonance analogue, see [HZel2]). 
To avoid dissipation of energy, we do not discuss the inverse resonance prob- 
lem in detail, but only mention some recent results closely related to the 
inverse spectral problems covered in this survey. We refer to [ZIO] for fur- 
ther discussion of inverse resonance problems for exterior domains and to 
[BJP, BP] for geometric scattering settings. 

The author would like to thank V. Baladi, G. Besson, Y. Colin de Verdiere, 
G. Courtois, R. Kuwabara, R. de La Llave, G. Lebeau, M. Rouleux, K. F. 
Siburg and M. Zworski for informative remarks on the contents of this survey. 
Of course, errors and omissions are the author's responsibility. 
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2 Expansions at t = 

The most classical spectral invariants are the heat invariants, namely the 
coefficients of the expansion at t = of the trace of the heat kernel. They 
are closely related to the coefficients of the trace of the wave group at t = 
0, although it should be noted that the powers t~2+"^ of the heat kernel 
expansion at t — are not singular if m is even and if m > | (where 
n — dimM). Hence, the heat kernel expansion contains more information 
than the singularity expansion of the wave trace at t — 0. 

2.1 Boundaryless case 

The earliest work in inverse spectral was based on calculations of heat in- 
variants in terms of curvature invariants, and the recovery of special metrics 
or domains from these curvature invariants [Ber, Pa, T, T2, Ku2, Ku3]. 
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We recall that the heat trace expansion in dimension n on a boundaryless 
manifold has the asymptotic expansion, 

oo 

i=o 

The coefficients aj are the heat invariants. We note that when n is odd, 
the powers of t are singular and hence the expansion may be viewed as a 
singularity expansion in which the terms become more regular. When n is 
even, the terms with —n/2 + j < are singular but the rest are smooth 
and hence are not residual. Rather one may view the expansion as a Taylor 
expansion at i = of r/2rre*^s. But just hke the singular terms, the 
coefficients are spectral invariants given by integrals of curvature invariants. 
The first four heat invariants in the boundaryless case are given by [T, T2] 

ao = Vol{M) = / dVolu 
ai = i J SdVolM 

(2.2) 

8 P'J T^kl Tfrs 8 T)rs jDjkl JD i 2 c| I?|2 

—^R^'^Rjj^Rijke — ^RjR'^Ri — ^S\Ric\^ + ^S^}dVolM- 

Here, S is the scalar curvature, Ric is the Ricci tensor and R is the Riemann 
tensor. In general the heat invariants are integrals of curvature polynomials 
of various weights in the metric. We refer to [Ber] for background. 

The heat invariants are complicated and it is difficult to detect mean- 
ingful patterns in the curvature polynomials. Nevertheless, they have been 
successfully used to obtain inverse spectral results. 

• Spheres: Tanno [T2] used 00,01,02,03 to prove that the round metric 
go on S"' for n < 6 is determined among all Riemannian manifolds 
by its spectrum, i.e. any isospectral metric g is necessarily isometric 
to go. He also used 03 [T2] to prove that canonical spheres are locally 
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spectrally determined (hence spectrally rigid) in all dimensions. Patodi 
proved that round spheres are determined by the spectra Spec^{M,g) 
and Spec^ (M, g) on zero and 1 forms. 

• Complex projective space: Let {M,g,J) be a compact Kahler man- 
ifold and let {CP"'{H), g^, Hq) be a complex n-dimensional projec- 
tive space with the Fubini-Study metric of constant holomorphic sec- 
tional curvature H. Tanno [T2] proves that if the complex dimension 
n < 6 and if Spec(M,5, J) = Spec(CP"(i/), ^o, ^o), then {M,g,J) 
is holomorphically isometric to {CP"'{H), g^, Jq). He also proves that 
{CP'^{H),gQ, Hq) is locally spectrally determined in all dimensions [T3]. 

• Flat manifolds: Patodi [Pa] and Tanno [T, T2] used the heat invariants 
to prove in dimension < 5 that if (M, g) is isospectral to a flat manifold, 
then it is flat. More precisely, they showed that if aj = for j > 1, 
and if n < 5 then {M,g) is fiat. The result is sharp, as Patodi (loc. 
cit.) showed that aj = for j > 1 for the product of a 3-dimensional 
sphere with a 3-dimensional space of constant negative curvature. In 
fact, Tanno showed that if a2 — = 0, then {M,g) is either E^/Ti, 
where Fi is some discontinuous group of translations of the Euclidean 
space E^ or (2) [S^C) x H^{-C)]/T2, where S%C)[H%-C)] is the 
3-spherc [hyperbolic 3-space] with constant curvature C > 0[—C < 0] 
and F2 is some discontinuous group of isometries of S^{C) x H^{—C). 
Kuwabara [Ku2, Ku3] used the invariants to prove that flat manifolds 
are locally spectrally determined, hence spectrally rigid. 

2.1.1 The boundary case 

When dfl ^ 0, the heat trace has the form 

00 

Tre'^^ - r"/' aj t^'^. (2.3) 
3=0 

The coefficients have been calculated for a variety of boundary conditions 
(see [BG] and its references). 

The formulae are simplest for plane domains, where the only invariant is 
the curvature k of the boundary. Using a nicely adapted calculus of pseudod- 
ifferential operators, L. Smith obtained the first five heat kernel coefficients 
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in the case of Dirichlet boundary conditions. They are given by 

Oo = area of fi, 



Oi = ~\f \ \ Q^\ (the length of the boundary), 



"2 — I /an 



< 



^3 = /sn 

^5 = /an '^^^^ ~ 2^ /an('^')^^*> 



(2.4) 



V 



Here, rfs is arclength and k/ is the derivative with respect to arclength. Cer- 
tain useful patterns in the heat coefficients were used by R. B. Melrose to 
prove a compactness result (later improved by Osgood-Philhps-Sarnak). We 
refer to [Me] for the details. In higher dimensions, one still has 



There exist a few inverse inverse spectral results using heat invariants: 

• Euclidean balls in all dimensions are spectrally determined among sim- 
ply connected bounded Euclidean domains by their Dirichlet or Neu- 
mann spectra. This follows from (2.5) and from the fact that isoperi- 
metric hypersurfaces in are spheres. 

• The exterior of the unit ball S3 C IR^ in dimension 3 is uniquely deter- 
mined among exterior domains of simply connected compact obstacles 
by its resonance poles. [HZ] . 

2.2 Domains and metrics with the same heat invari- 



It was soon realized that heat invariants are insufficient to determine smooth 
metrics or domains. This is due to the fact that they are integrals of local 



(2.5) 



ants 
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invariants of the metrics. Pairs of non-isometric metrics with the same heat 
invariants can be obtained by putting two isometric bumps. The bumped 
spheres will not be isometric if the distances between the bumps are different, 
but the heat invariants will be the same. There are many variations on this 
well-known example. But heat invariants might be quite useful for analytic 
metrics and domains, and have also been used in compactness results. 



3 Dynamics and dynamical inverse problems 

We now turn to the more dynamical theory of the wave group in inverse 
spectral theory. The trace of the wave group expresses spectral invariants 
in terms of the dynamics of the geodesic flow, and often 'reduces' inverse 
Laplace spectral problems to inverse problems in dynamics. We therefore 
begin by recalling the relevant dynamical notions and inverse problems. 



3.1 Geodesic flow on boundaryless manifolds 

We denote by {T*M, Ylj dxj A d^j) the cotangent bundle of M equipped 
with its natural symplectic form. Given a metric g, we define the metric 
Hamiltonian in a standard notation by 



n+l 



E^''(^)^^^^- (3-1) 

ij=l 



and define the energy surface to be the unit sphere bundle S*M — {{x,^)\ : 
l^lg = 1}. The spectral theorists geodesic flow is the Hamiltonian flow 

= exptEn : T*M\ T*M\0, Eh = the Hamiltonian vector field of H. 

(3.2) 

It is homogeneous of degree 1 with respect to the dilation (x, ^) — > (x, r^), r > 
0, so nothing is lost by restricting G* to S*M. We also denote its generator 
by S, the Hamiltonian and metric being understood. 

The periodic orbits of of period T of the geodesic fiow are the fixed 
points of on S*M. Equivalently they are the critical points of the length 
functional on the free loop space of M of length T. 
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3.2 Billard flow and billiard map on domains with 
boundary 

To define the geodesic or billiard flow G* on a domain Q with boundary d^, 
we need to specify what happens when a geodesic intersects the boundary. 
The definition is dictated by the propagation of singularities theorem for 
solutions of the wave equation (due to R. B. Melrose and J. Sjostrand [MS]), 
and is therefore not purely geometric or dynamical. A billiard trajectory is 
defined to be the path along which a singularity of a solution of the wave 
equation propagates. We give a quick, informal review of the fiow; for more 
details, the reader might consult [MS], [PS]. For simplicity we assume that 
there are no points of infinite order tangency. 

We denote by S*fl the unit tangent vectors to the interior of Q and by 
S*^dD, the manifold with boundary of inward pointing unit tangent vectors 
to Q with footpoints in dQ. The boundary consists of unit vectors tangent 
to dft. The biUiard fiow G* is a flow on S*Q U S*j^dQ, deflned as follows: 
When an interior geodesic of Q intersects the boundary dQ transver sally, 
it is reflected by the usual Snell law of equal angles. Such trajectories are 
called (transveral) reflecting rays. The complications occur when a geodesic 
intersects the boundary tangentially in S*dQ. 

Convex domains arc simpler than non-convex domains, since interior rays 
cannot intersect the boundary tangentially. Dynamical studies of billiards 
(see e.g. [MF]) often restrict to convex domains. It should be noted that 
geodesies of with the induced metric are important billiard trajectories. 
They are limits of 'creeping rays', i.e. rays with many small links (interior 
segments) which stay close to dQ. In particular, the boundary of a convex 
plane domain is a closed billiard trajectory. In higher dimensions, closed 
geodesies on the boundary which are limits of interior creeping rays are closed 
biUiard trajectories. (Are all closed geodesies on the boundary limits in this 
sense?) 

On a non-convex domain, a general billiard trajectory is divided into seg- 
ments which are either geodesic segments in the interior or geodesic segments 
of the boundary. Geodesic segments of the boundary only occur where the 
boundary is convex. Intuitively, the boundary segments are limits of creeping 
rays along the convex parts of the boundary. Trajectories enter and exist the 
boundary at inflection points. In particular, the boundary of a non-convex 
plane domains is not a closed billiard trajectory. If a trajectory intersects 
the boundary tangentially at a non-inflection point of a non-convex domain. 
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it goes straight past the point of intersection into the interior. 

When {Q. g) is non-trapping (i.e. if there is no geodesic ray which remains 
forever in the interior), the set S*^dQ behaves hke a global cross section to 
the billiard flow. It is then natural to reduce the dimension by defining the 
billiard ball map p : B*{dn) B*{dn), where B*{dn) is the ball bundle 
of the boundary. We first identify S*^dO. ~ B*{dO.) by adding to a tangent 
(co)vector t] G B*dQ of length < 1 a multiple cu^ of the inward point unit 
normal Ug to form a covector in S^^Q. The image f3{q,v) is then defined to 
be the tangential part of the first intersection of G*{q, rj + cuq) with dfl. The 
billiard map is symplectic with respect to the natural symplectic form on 
B*{dQ). 

An equivalent description of the billiard map of a plane domain is as 
follows. Let q G dQ and let G (0,7r). The point {q,4>) corresponds to an 
inward pointing unit vector making an angle (p with the tangent line, with 
= corresponding to a fixed orientation (say counter-clockwise). The 
billiard map is then P{q, (f) — [q' , (/)') where {q\ 0') are the parameters of the 
refiected ray at the first point of intersection with the boundary. The map f3 
is then area preserving with respect to sin 4>ds A d4> (see e.g. [MF]). 

3.3 Closed orbits and their Poincare maps 

Closed orbits (or periodic orbits) 7 of fiows are orbits of points {x, ^) E T*M 
satisfying G-'^{x,^) = (x,^) for some T 7^ (the period). They project to 
closed geodesies on the Riemannian manifold or domain. 

We recall the definition of the nonlinear Poincare map 7^- in S*M one 
forms a symplectic transversal to 7 at some point mo- One then defines 
the first return map, or nonlinear Poincare map, 

•^7(0 • 'S'-y > S.y 

by setting = G'^^'^\C), where T{Q is the first return time of the tra- 

jectory to S^. This map is well-defined and symplectic from a small neigh- 
borhood of 7(0) = mo to a larger neighborhood. By definition, the linear 
Poincare map is its derivative, Py = d'?^{m.o). 

Closed geodesies are classified by the spectral properties of the symplectic 
linear map P^. Its eigenvalues come in 4-tuples A,A, A~^,A~^. A closed 
geodesic 7 is called: 

• non- degenerate if det(/ — P^) ^ 0; 
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• elliptic if all of its eigenvalues are of modulus one and not equal to ±1, 
in which case they come in complex conjugate pairs e*^°^ . 

• hyperbolic if all of its eigenvalues are real, in which case they come in 
inverse pairs XjXj^ 

• loxodromic or complex hyperbolic in the case where the 4-tuple consists 
of distinct eigenvlaues as above. 

There are other possibilities (parabolic) in the degenerate case. In the 
case of Euclidean domains, or more generally domains where there is a unique 
geodesic between each pair of boundary points, one can specify a billiard 

trajectory by its successive points of contact qo, qi, q2, ■ ■ ■ with the boundary. 
The n-link periodic reflecting rays are the trajectories where g„ = go for some 
n > 1. The point qo, ... ,qn is then a critical point of the length functional 

n-l 

^{qo,---,qn) = ^\qi+i- qi\ 

i=0 

on {dn)"". 

Among the periodic orbits, bouncing ball orbits often have special ap- 
plications in inverse spectral theory. By a bouncing ball orbit 7 one means 
a periodic 2-link reflecting ray, i.e. go — 0.2- Convex domains always have 
at least two bouncing ball orbits, one of which is its diameter. There exist 
non-convex domains without any bouncing ball orbits. For geometric aspects 
of bouncing ball orbits, we refer to [Gh]. 

The projection to Q consists of a segment gogl which is orthogonal to the 
boundary at both endpoints, i.e. gogl is an extremal diameter. The period is 
of course twice the length of the segment, which we denote by L. We write 
go = A,qi = B. We orient the domain so that the links of 7 are vertical and 
so that the midpoint is at the origin of R^. The top, resp. bottom, of the 
boundary is then the graph of a function y — f+{x) resp. y — f-{x) over 
the X-axis. It is clear that the wave invariants depend only on the Taylor 
coefficients of f± at A, B. 

We denote by Ra, resp. Rb the radius of curvature of the boundary at 
the endpoints A, resp. B of an extremal diameter. The bouncing ball orbit 
is elliptic if L < min{i?^, Rb} or if max{i?^, Rb} < L < Ra -\- Rb, and is 
hyperbohc if L > Ra + Rb or if mm{RA,RB} < L < miix{RA,RB} (see 
[KT]). If AB is a local minimum diameter than L < Ra + Rb while if it is 
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a local maximum diameter then L > Ra + Rb- So a (non-degenerate) local 
maximum diameter must be hyperbolic; a local minimum diameter is elliptic 
if it satisfies the additional inequalities above. 

When 7 is elliptic, the eigenvalues of are of the form {e^*"} while in 
the hyperbolic case they are of the form {e^^}. The explicit formulae for 
them are: 

cosq;/2 = ^ — ;^)(1 ~ '^~) (elliptic case), 

" (3.3) 
cosh A/2 = ^(1 - ;^)(1 - (hyperbolic case). 

We note that /±(0) = ^ = k±, where k± denotes the curvature at (^4 = 
-,5 = +)- 

3.4 The length spectrum and the marked length spec- 
trum 

The length spectrum of a boundaryless manifold (M, g) is the discrete set 

Lsp{M, g) = {L^, < L^, < • • • } (3.4) 

of lengths of closed geodesies 7^-. In the boundary case, the length spectrum 
Lsp{Q) is the set of lengths of closed billiard trajectories in the sense of §3.2 
and is no longer discrete, but rather has points of accumulation at lengths 
of trajectories which have intervals along the boundary. In the case of con- 
vex plane domains, e.g., the length spectrum is the union of the lengths of 
periodic reflecting rays and multiples of \dQ\. According to the standard ter- 
minology, Lsp{M, g) is the set of distinct lengths, not including multiplicities, 
and one refers to the the length spectrum repeated according to multiplicity 
as the extended length spectrum. 

In the notation for Lsp{M, g) we wrote L^^ as if the closed geodesies of 
this length were isolated. But in many examples (e.g. spheres or flat tori), 
the geodesies come in families, and the associated length T is the common 
length of closed geodesies in the family. In place of closed geodesies, one has 
components of the fixed point sets of G^^ at this time. The fixed point sets 
could be quite messy, so it is also common to assume that they are clean, 
i.e. that the fixed point sets are manifolds, and that their tangent spaces 
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are fixed point sets of dG^ . It is equivalent tfiat the lengtli functional is 
Bott-Morse on the free loop space. 

The set of lengths is unformatted in the sense that one does not know 
which lengths in the list correspond to which closed geodesies. A formatted 
notion in which lengths are assigned to topogically distinct types of closed 
geodesies is the marked length spectrum MLg. On a manifold without bound- 
ary, it assigns to each free homotoply class of closed loops the length of the 
shortest closed geodesic in its class. On a convex plane domain with bound- 
ary, topologically distinct closed geodesies correspond to different rotation 

numbers — = — winding number — rpj^^ marked length spectrum then as- 

" number qt reflections , , , r , 

sociates to each rational rotation number the maximal length of the closed 

geodesies having n reflection points and winding number m. We refer to §8.4 

(based on [SI, S2]) for further discussion and applications to inverse spectral 

theory. 



3.5 Birhoff normal forms 

Birkhoff normal forms are approximations to Hamiltonians (or symplectic 
maps) near equilibria by completely integrable Hamiltonians (or symplectic 
maps). We now briefly consider the simplest kind of equilibrium occurring 
for geodesic flows, a closed geodesic 7. 

Let us first consider the Birkhoff normal form of the metric Hamiltonian 
(3.1) near a non- degenerate elliptic closed geodesic 7. To put H into normal 
form is to conjugate it (approximately) to a function of locally defined action 
variables (cr, /i, . . . , /„) on the model space T*{S^ x M"), where a is the 
momentum coordinate in T*S^ and Ij are transversal action variables which 
depend on the type (elliptic, hyperbolic, loxodromic) of 7. More precisely, 
the normal form algorithm defines a sequence of canonical transformations 
Xm at 7 which conjugate H to the normal forms 

XmH = cr+-rZ^ oijlj + + ••• + m ^m+i 

L ^ — ' (7 (7 

(3.5) 

where is homogeneous of order k+1 in Ji. . . . , /„, and where 0\[j^^ is the 
space of germs of functions homogeneous of degree 1 which vanish to order 
M + 1 along 7. Note that all the terms in (3.5) are homogenous of degree 1 
in (cr, /i, /„), and that the order of vanishing at |/| = equals one plus the 
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order of decay in a. The coefficients of the monomials in the Pj{Ii, ■ ■ ■ , In) 
are known as the classical Birkhoff normal form invariants. 

The algorithm for putting a symplectic map (or Hamiltonian) into Birkhoff 
normal form around a fixed point (or equilibrium point) can be found in 
many places (see e.g. [AKN]). For geodesic flows, we need Birkhoff normal 
forms for homogeneous Hamiltonians around periodic orbits. An algorithm 
for putting a metric Hamiltonian into Birkhoff normal form around a closed 
geodesic is described in [G] and also in the appendix to [Z3], among other 
places. To carry out the algorithm to infinite order, one needs to assume that 
no eigenvalues of 7 are roots of unity, which is of course a stronger condition 
than non-degeneracy. 

There are also a Birkhoff normal forms for the geodesic flow and for the 
Poincare map. The normal form is simpler for since we have eliminated 
the directions along 7 and therefore may express just in terms of local 
action-angle variables (/, 0). The normal form is as follows: 

7{I,cf>)^{I,cf> + WiGM{I)), mod Oi,+i, (3.6) 

where Gm{I) is a polynomial of degree M in the / variables. Thus, to 
order M + 1, the Poincare map leaves invariant the level sets of the actions 
(ellipses, hyperbolas etc. according to the type of 7) and 'rotates' the angle 
along them. 

The well-known question arises whether the full (infinite series) Birkhoff 
normal form converges and whether the formal symplectic map conjugating 
the Hamiltonian to its normal form converges. The latter is sometimes called 
the Birkhoff transformation. Clearly, if the Birkhoff normal form diverges, 
then so must the Birkhoff transformation. According to a recent article of 
Perez-Marco [PM], there are no known examples of analytic Hamiltonians 
having divergent Birkhoff normal forms, and we refer to that article for fur- 
ther background and history on the problem. If both the normal form and 
the transformation converge, the Hamiltonian must be integrable in a neigh- 
borhood of the orbit. The Birkhoff normal form is then the expression of 
the Hamiltonian in local action-angle variables. In the generic non-resonant 
case, it was proved by H. Ito that the Birkhoff transformation converges [II] 
(see also [12] for the resonant case). 

For non-integrablc systems, the Birkhoff normal form (3.5) and the con- 
jugating map an around elliptic orbit 7 are only approximations in small 
neighborhoods of 7, which shrink at M increases. The hyperbolic case is 
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simpler. In the case of hyperbolic orbits of analytic symplectic maps in two 
degrees of freedom, it was proved by J. Moser [Mo] that the Birkhoff nor- 
mal form and transformation do converge. Convergence is more complicated 
in higher dimensional, and we refer to Banyaga-de La Llave- Wayne [BLW], 
Perez-Marco [PM] and to Rouleux [R] for recent results. In Rouleux, results 
of [BLW] are used to prove the existence of a local smooth symplectic conju- 
gacy K*H = q{I) of the metric Hamiltonian H near a hyperbolic fixed point 
(or orbits) to a smooth normal form obtained by Borel summation of 

the formal Birkhoff normal form of H. (In the elliptic case, one only has a 
conjugacy k,*H = q{I) +r up to a remainder r which vanishes to infinite order 
at the orbit.) It follows that two metric Hamiltonians with the same Birkhoff 
normal form at respective hyperbolic closed geodesies have locally symplec- 
tically equivalent geodesic flows near those orbits. (The author thanks M. 
Rouleux for corroborating this point). 

3.6 Livsic cohomology 

The cohomology problem asks whether a function (cocycle) F G C°°(S'*M) 
satisfying Fds = for every closed geodesic of the metric g is necessarily 
a co-boundary, F — S(/) where S is the generator of the geodesic flow G* 
and / is a function with some degree of regularity (see [LI] and other works 
of de la Llave and others for regularity results). 

It is relevant to the inverse length spectral problem for the following 
reason, first observed by Guillemin-Kazhdan [GK]: Under a deformation of 
a metric g = go preserving the extended Lsp{M, g) (including multiplicities), 
one has 



When the cohomology is trivial, one can therefore write g = S(/) for some 
/ with the given regularity. Then one can study the harmonic analysis on 
S*M to see whether the symmetric 2-tensor g can be expressed as S(/). The 
answer is no for negatively curved surfaces, and that gave the rigidity result 
of [GK]. Since then their result has been improved (see [CS] for the most 
general result on rigidity for negatively curved manifolds), but the general 
strategy retains some potential for other kinds of metrics. 

It is known that the Livsic cohomology problem is always solvable for an 
Anosov flow on a closed manifold. We refer to [LI, Wa] for recent results. One 




7 



(3.7) 
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might ask for the analogous result for hyperbolic billiard flows on bounded 
domains. It is difficult to formulate the analogous regularity result since the 
dynamics are not smooth. For instance, the billiards are not even defined at 
the corners of a domain with piecewise smooth boundary and with concave 
boundary faces. One studies billiards on the domain by puncturing out the 
measure zero set of orbits which ever run into the corners. The resulting phase 
space is then not a closed manifold and there would be complicated issues 
about regularity of solutions as one approached the punctured set. Perhaps 
the simplest setting would be that of Sinai billiards, i.e. the exterior of a 
convex obstacle in a compact manifold (e.g. the exterior of a disc in a two- 
dimensional torus). The billiards are hyperbolic and there are no glancing 
orbits or corner orbits to puncture out. But we are not aware of any work 
on the Livsic equation in this setting. 

The Livsic equation can be studied for for non-hyperbolic flows, but there 
do not seem to exist many studies of it other than manifolds all of whose 
geodesies are closed (in which case a simple Fourier analysis suffices). One 
study in the near-integrable case is the article [L2] of de la Llave. One does 
not expect the cohomology to be trivial in general settings, but the results 
on this equation might have interesting implications for the length spectral 
deformation problem. As will be explained in §4.3, in certain cases such 
as integrable systems, the relevant homological equation is not the Livsic 
equation but an analogue where one integrates over non-degenerate critical 
manifolds of geodesies. 

3.7 Dynamical inverse problems 

As mentioned above, inverse spectral theory often proceeds by showing that 
certain dynamical invariants are invariants of the Laplace spectrum. This 
gives rise to inverse dynamical problems. In this section, we briefly survey 
some of the problems and results. A recent survey of rigidity and conjugacy 
problems has been written by C. Croke [C3]. We are grateful to V. Baladi, 
G. Besson, G. Courtois, C. Gordon, R. de La Llave, and K. F. Siburg for 
advice on this section. 

One of the important problems in dynamics is the following inverse spec- 
tral problem in dynamics. 

Problem 3.1 For which {M,g) does the marked length spectrum of{M,g) 
determine {M,g) up to isometry? 
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In the positive direction, V. Bangert [B] proved that the marked length 
spectrum of a flat two-torus determines the flat metric up to isomctry. J. P. 
Otal [0] and C. Croke [C] (see also [CFF, C3]) independently proved that the 
marked length spectrum determines surfaces of negative curvature; for more 
general results sec [C3]. U. Hamenstadt [H] proved that a locally symmetric 
manifold is determined by its marked length spectrum. An example due 
to F. Bonahon shows that metric structures more general than Riemannian 
metrics are not always determined by their marked length spectra. 

A related problem concerns the conjugacy rigidity of Riemannian mani- 
folds. The geodesic flows G* of Riemannian manifolds (Mj, gj) are called C*^ 
conjugate if there exists a time-preserving C*^ conjugacy between them, i.e. 
a C'^ homeomorphism x '■ ^gi^i ~^ ^g2^^ satisfying 

XoG\ox-' = Gl. (3.8) 

If X extends to a homogeneous symplectic diffeomorphism of the cotangent 
bundles, the flows are called symplectically conjugate. Two billiard maps are 
conjugate if there exists a symplectic diffeomorphism x '■ B*D,i — > B*D,2 such 
that 

XO/9iox"' = /32. (3.9) 

Problem 3.2 When does existence of a conjugacy between two geodesic 
flows (for a given k) imply isometry of the metrics? Does symplectic equiv- 
alence of billiard maps of convex domains imply their isometry? 

In some cases, equality of marked length spectra implies C° conjugacy of 
geodesic flows. It was proved independently by J. P. Otal [0] and by C. Croke 
[C] that negatively curved surfaces with the same marked length spectrum 
have C°- conjugate geodesic flows, and that such surfaces must be isometric. 
The corresponding statement in higher dimensions still appears to be open 
in general. It is stated in [C3] that compact flat manifolds are C°° conjugacy 
rigid. In [H2], Hammenstadt constructed a time-preserving conjugacy 
between negatively curved manifolds with the same marked length spectrum. 
When this conjugacy is the work of Besson, Courtois and Gallot [BCG] 
proves that the manifolds must be isometric. The method of [H] avoids 
this regularity issue. The conjugacy problem for certain nilmanifolds was 
studied by C. Gordon, D. Schueth and Y. Mao in [GorM, GMS]. They proved 
for special classes of 2-step nilmanifolds that C° conjugacy of the geodesic 
flow implies isometry. These manifolds occur in a non-trivial isospectral 
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deformation and have the same marked length spectrum, so these examples 
also show that equality of marked length spectra does not necessarily imply 
C° conjugacy of the geodesic flows. 

There is at least one known example where symplectic conjugacy does 
not imply isometry, namely for surfaces all of whose geodesies are closed. 
The geodesic flow of such a ZoU surface was shown by Weinstein to be sym- 
plectically equivalent to that of the standard 2-sphere. Aside from related 
examples, all of the other results known to the author (see [C3]) are rigid- 
ity results showing that conjugacy implies isometry. The relations between 
marked length spectral equivalence, conjugacy of geodesic flows and isometry 
do not seem to have been studied in other settings. 

One could pose local versions of the symplectic conjugacy problem, where 
there should be much less rigidity. A local version (on the level of germs 
around closed orbits) of symplectic conjugacy is symplectic conjugacy be- 
tween Poincare maps: 

x-s,,G s;^M, ^ s^, c s;^M2, x'S'^.x'^ = ^72- (3.10) 

One could ask: 

Problem 3.3 When does local symplectic conjugacy of Poincare maps at a 
closed geodesic (or local symplectic conjugacy of geodesic flows) imply local 
isometry? 

As a special case of the local equivalence problem, suppose that (M, g) 
is a hyperbolic manifold, and let 7 be a closed geodesic. Let g' be a second 
real analytic metric on M. Suppose that there exists a closed geodesic for g' 
for which the Poincare maps CP^ are symplectically conjugate. Must g' be a 
hyperbolic metric? 

Even less rigid is the Birkhoff normal form, since in general there only 
exists a formal (power series) canonical transformation conjugating the germs 
of the metrics. 

Problem 3.4 To what extent can the germ of a metric g at^ be determined 
by its Birkhoff normal form at j? 

The sceptic might suspect that when the formal conjugacy does not con- 
verge, the Birkhoff normal form gives little beyond a list of numerical in- 
variants, similar to and no simpler than the wave invariants. When it does 
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converge, there might be too httle rigidity to determine very much about the 
metric. 

For instance, it was observed by Y. CoUn de Verdiere [CdV] the Birkhoff 
normal form of the bilhard map around a bouncing baU orbit of a domain 
with the symmetries of an eUipse determines all of the Taylor coefficients of 
the domain at the endpoints of the orbit. Hence a real analytic domain with 
two symmetries is determined by the Birkhoff normal form of the Poincare 
map at a bouncing ball orbit. 

If the domain has only one symmetry, then the calculation of [CdV] shows 
that the Birkhoff normal form at a bouncing ball orbit does not contain 
enough information to determine the Taylor coefficients at both the top and 
bottom of the domain. This simple, 2-dimensional example shows how weak 
an invariant the Birkhoff normal form is. 



3.7.1 Zeta functions of geodesic flows 

We add a final inverse problem which to our knowledge has not been studied 
before. Let us consider two zeta functions associated to the geodesic flow. 
The first is the Ruelle zeta function defined by 

oo —skL{-y) 

Here, we are assuming that the metric is bumpy so that the fixed point sets 
of the geodesic flow consist of isolated, non-degenerate closed orbits. We 
denote by CP = {7} the set of primitive periodic orbits of the geodesic flow 
and by Py the linear Poincare map of 7. As will be seen below, it is a spectral 
invariant as long as the (extended) length spectrum is simple (multiplicity 
free) or if the metric has no conjugate points. 

The zeta function log Z(^M,g){s) arises as the Laplace transform of the 
so-called flat trace of the unitary Koopman operator 

Wt : L'{S*M) ^ L'{S*M), WJ{x,0 ■= f{G\x,0)- 

This operator does not possess a distribution trace in the usual sense, but 
the integral of its diagonal part Jg,j^jW{t,uj,uj)dV{uj) is well-defined in the 
sense of pushing forward and pulling back distributions. In this sense, it was 
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observed by V. Guillemin [G2] that 



TrWit) = C.VoliM,g)S^-\t) + EE ^^T^' (3-12) 

The Laplace transform may be viewed as the flat trace of the resolvent of 5, 
so that 

log Z(M,5)(s) = Tr(S-s)-\ (3.13) 

In a formal sense, the Ruelle zeta function is a trace of a spectral function 
of the geodesic flow and it is natural to wonder how much of the spectrum of 
Wt can be determined from the traces. In particular, crgodicity, weak mixing 
and mixing are spectral properties of geodesic flows, i.e. can be read off from 
the spectrum of Wf. For instance, ergodicity is equivalent to the statement 
that the multiplicity of the eigenvalue 1 for Wt equals one. 

But the trace is not a distribution trace on in the sense that if p e C^, 
f^p(t)Wtdt is usually not a trace class operator on L'^{S*M). Hence, it is 
far from clear that one can read off spectral properties of Wt from these 
traces. At this time of writing, almost nothing seems to be known about 
zeta functions except in the case of hyperbolic flows, in which case one knows 
many properties of the spectrum of the geodesic flow, e.g. that it is mixing. 
We refer to papers of V. Baladi (see e.g. [Ba] for her survey) for background 
on spectral interpretations of zeta functions and to C. Liverani [BKL] et. al. 
for some possible future extensions. C. Deninger and others (sec e.g. [DS]) 
have recently studied the flat trace in non-hyperbolic dynamical settings. 

It might be interesting to develop the theory of zeta functions for broader 
classes of flows. Indeed we are interested in the inverse problem of determin- 
ing the dynamical properties from the trace. 

A related zeta function is deflned for 3?s large by 

It is a Laplace spectral invariant of metrics without conjugate points, even if 
(as with hyperbolic quotients) the extended length spectrum has multiplicity. 

Problem 3.5 What are the analytical properties of these zeta functions 
when the geodesic flow has some given dynamical signature. Can one de- 
termine apriori from Z(^M,g){s) whether the geodesic flow is ergodic, weak 
mixing or mixing? 
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The first step is to see how far the zeta function can be analytically 
continued, and whether one can determine the multiplicity of the pole at 
s = 0. We then ask if the multiplicity has a spectral interpretation as the 
dimension of the invariant L^-functions under the flow. 

We observe that the zeta functions C(M,g) (s) are Laplace spectral invariant 
for metrics with multiplicity free length spectra as well as for those without 
conjugate points. Hence this dynamical inverse result would immediately 
imply that one can determine spectral properties of the geodesic flow from 
the Laplace spectrum. 

4 Wave invariants 

We now return to the Laplace spectrum and introduce the wave invariants, 
which are a more discriminating set of spectral invariants than heat invari- 
ants. They arise from the trace of the wave group U{t) — e**^ of [M,g). 
One forms the (distribution) trace 

TrU{t) = (4-1) 

It is a tempered distribution on IR. We denote its singular support (the 
complement of the set where it is a smooth function) by Sing Supp TrU (t) . 

The first result on the wave trace is the Poisson relation on a manifold 
without boundary. 

Sing SupprrC/(i) C Lsp(M,g), (4.2) 

proved by Y. Colin de Verdiere [CdV2, CdV3], Chazarain [Ch2], and Duistermaat- 
Guillemin [DG] (following non-rigorous work of Balian-Bloch [BB2] and Gutzwiller 
[Gutz]). The generalization to manifolds with boundary was proved by 
Anderson-Melrose [AM] and Guillemin-Melrose [GM]. As above, we denote 
the length of a closed geodesic 7 by L^. For each L = G Lsp{M,g) there 
are at least two closed geodesies of that length, namely 7 and 7"^ (its time 
reversal). The singularities due to these lengths are identical so one often 
considers the even part of TrU{t) i.e. TrE{t) where E{t) — cos{ty/A.). 

We emphasize that (4.2) is only known to be a containment relation. As 
will be seen below, cancellations could take place if a length L G Lsp{M, g) 
is multiple, so that TrU{t) might be smooth at L G Lsp{M,g). 
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Problem 4.1 Is Lsp{M,g) a spectral invariant? Are there any examples 
where TrU{t) is C°° att^L^^ Lsp{M, g) ? 

No example of this seems to be known. Y. Colin de Verdiere has pointed 
out that it is even unknown whether TrU (t) could be smooth on all of M\{0}. 
Moreover, Lsp{M, g) docs not include information about multiplicities of 
lengths. Sunada-type isospectral pairs always have multiple length spectra, 
and for many (presumably, generic) examples, the length spectra have differ- 
ent multiplicities. The reason is that the geodesies on the base manifold of a 
Sunada quadruple can split in different ways in the covers. Other examples 
of isospectral pairs (both in the sense of eigenvalue and length spectra) with 
different multiplicities of lengths have been constructed by R. Cornet, C. 
Cordon and others. 

The indeterminacy of multiplicities raises a natural problem. Let us recall 
that the topological entropy of the geodesic flow of (M, g) is the exponential 
growth rate of the length spectrum: 

htop ■■= liminf log#{7 : < L}. (4.3) 

L— »oo 

Here, lengths are counted with multipicity. In the case where geodesies come 
in families, we count components of the fixed point sets. 

Problem 4.2 Is htop cl spectral invariant? 

The question is only non-trivial when multiplicities in the length spec- 
trum grow as fast as the length spectrum, as occurs for compact hyperbolic 
manifolds in dim > 3 and for arithmetic hyperbolic quotients in dimension 
2 . As observed by Besson-Courtois-Callot [BCC], an affirmative answer 
would show that that hyperbolic metrics are spectrally determined among 
other negatively curved metrics, since they are the unique minimizers of htop- 
A closely related problem is: 

Problem 4.3 (see [BCG]) Suppose that M with dimM > 3 possesses a 
hyperbolic metric, and let M_ denote the class of negatively curved metrics 
on M. Is the hyperbolic metric g determined by its spectrum among metrics 
in M_ ? I.e. can there exist another non-isometric metric in this class which 
is isospectral to the hyperbolic metric? 
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4.0.2 SingulEir support versus analytic singuleir support 

While discussing the Poisson relation, we pose the following question con- 
cerning the analytic Poisson relation: 

Problem 4.4 Can one tell from S'pec(A) if a metric (or the underlying 
manifold or domain) is real analytic? 

The idea is to calculate the analytic wave front set of TrU{t), i.e. 

The analytic wave front set is the complement of the set where the trace is 
real analytic. When (M, g) is a (real analytic) Riemannian manifold, the 
analytic wave front set WFaTrU{t) is also the set Lsp{M, g), that is, TrU{t) 
is a real analytic function outside of this discrete set. If (M, g) is a C°° but 
not a C"^ Riemannian manifold, it is plausible that WFaTrU(t) could contain 
an interval or be all of M. Thus simply the discreteness of WFaTrU {t) would 
say that (M, g) is real analytic. 

4.1 Singularity expansions 

Much more is true than the Poisson relation: TrU{t) has a singularity ex- 
pansion at each L e Lsp{M,g): 

TrU{t) = eo{t) + Y.LeLsv(M,g) ^L{t) mod (4.4) 

where cq, arc Lagrangean distributions with singularities at just one point, 
i.e. singsuppeo = {0}, singsuppcL = {L}. When the length functional on the 
loopspace of M is a Bott-Morse functional, the terms have complete asymp- 
totic expansions. In the Morse case (i.e. bumpy metrics), the expansions 
take the form 

eo{t) = ao,-n{t + iO)-" + ao,-n+i(t + iO)-"+^ + • • • (4.5) 
eiit) = aL,-i{t-L + iO)-' + aL,olog{t-{L + iQ)) 

(4.6) 

+ aL,i{t-L + iO)\og{t-{L + iO)) + --- , 
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where • • • refers to homogeneous terms of ever higher integral degrees ([DG]). 

The wave coefficients ao,fe at t = are essentially the same as the singular 
heat coefficients, hence are given by integrals over M of /^^ PjiR, Vi?, ...)dvol 
of homogeneous curvature polynomials. The wave invariants for t ^ have 
the form: 



where a^j involves on the germ of the metric along 7. Here, {7} runs over 
the set of closed geodesies, and where L^, L*, m^, resp. are the length, 
primitive length, Maslov index and linear Poincare map of 7. (The primitive 
length of a closed orbit is the least non-zero period of the orbit, i.e the length 
once around). For instance, the principal wave invariant at t — L in the case 
of a non-degenerate closed geodesic is given by 



The same formula for the leading singularity is valid for periodic reflecting 
rays of compact smooth Riemannian domains with boundary and with Neu- 
mann boundary conditions, while in the Dirichlet case the numerator must be 
multiplied by (—1)^ where r is the number of reflection points (see [GM, PS]). 

The wave invariants for t 7^ are both less global and more global than 
the heat invariants. First, they are more global in that they are not integrals 
of local invariants, but involve the semi-global first return map !P^. One 
could imagine different local geometries producing the same first return map. 
Second, they are less global because they arc determined by the germ of the 
metric at 7 and are unchanged if the metric is changed outside 7. 

Thus, associated to any closed geodesic 7 of (M, g) is the sequence {aYj} 
of wave invariants of 7 and of its iterates Y- These invariants depend only 
on the germ of the metric at 7. The principal question of this survey may 
be stated as follows: 

Problem 4.5 How much of the local geometry of the metric g at j is con- 
tained in the wave invariants {a^j} ? Can the germ of the metric g at ^ be 
determined from the wave invariants? At least, can the symplectic equiva- 
lence class of its germ be determined? 

As will be discussed in the next section §6, one of the principal results on 
this problem is the theorem (due to V. Guillemin [G, G2], with some input 




a. 



(4.7) 



'y:L-y=L 




(4.8) 



33 



and generalizations by the author [Z3, Z4]) that the classical Birkhoff normal 
form of the metric (or the Poincare 7 is determined by the wave 

trace invariants. On a global level: 

Problem 4.6 How much of the global geometry (M, g) is contained in the 
entire set of wave invariants {a^ j}? 

As these questions suggest, one may divide the potential use of wave 
invariants into two classes: (i) those which use all of the closed geodesies, 
and (ii) those which involve one or a few closed geodesies. Obviously, (i) 
is more powerful if one can combine information from all of the geodesies, 
since it adds precisely the global feature which is lost by studying just one 
closed geodesic. But it seems very difficult in general to combine information 
coming from different geodesies. 

To the author's knowledge, the global problem of combining wave in- 
variants of all closed geodesies has only been led to successful results on 
isospectral deformations. We will briefly survey the methods and results in 
the next section §4.3. Otherwise, the main results use only one or two closed 
geodesies, and this cannot possibly succeed unless the metrics or domains 
arc assumed real analytic. We will survey the results in the analytic case in 
§7.5. 

Another dichotomy in the use of wave invariants is whether one uses only 
'principal term' information at each geodesic, i.e. the invariant 07,-1 (4.8), or 
whether one uses all of the terms. Before the latter is possible, one needs to 
calculate the lower order terms to some degree. We will discuss the possible 
calculations in detail in later sections. 

4.2 Inverse spectral results using wave invariants 

We now consider results which use wave invariants to specify a class of metrics 
or domains which are spectrally determined in the class or which admit no 
isospectral deformations. We will discuss the proofs of some of the results in 
later sections to illustrate the methods. 

The positive results based on wave invariant analysis are as follows. 

• Negatively curved compact manifolds are spectrally rigid [GK, C]. 

• Simple real analytic surfaces of revolution of 'simple type' (with one 
critical distance from the axis) are spectrally determined within the 
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class of such surfaces [Z6]. Any other surface which is isospectral to 
a simple surface of revolution must be C°-integrable [Z7, S2]. Smooth 
surfaces of revolution with a mirror symmetry thru the x — y plane are 
spectrally determined among metrics of this kind [BH]. 

• Simply connected analytic plane domains with two symmetry axes (i.e. 
with the symmetries of an ellipse) and with a bouncing ball orbit of 
fixed length L are spectrally determined within this class ([Zl, ISZ, 
GM2], see also [CdV] for an earher result proving spectral rigidity of 
domains in this class). A closely related result is that convex analytic 
domains with two symmetry axes are spectrally determined within this 
class [Zl]. The shortest orbit is necessarily a bouncing ball orbit and 
of course its length is a spectral invariant [Gh]. 

• Simply connected analytic plane domains with one symmetry, and with 
a bouncing ball orbit whose orientation of a fixed length L which is 
reversed by the symmetry, are spectrally determined within this class 
[Z5]. (This implies the preceding result, but we state it separately since 
it is a new result based on different methods). 

• There is a spectrally determined class of convex plane domains (el- 
lipses?) for which each element is spectrally determined among all 
convex plane domains [MM]. (Other isolated spectrally determined 
examples have recently been given in [W]). 

• The mean minimal action of a convex billiard table is invariant under 
isospectral deformations [SI] (see §8.4 for further discussion). 

• The exterior of a two-component analytic obstacle with two symmetries 
around a bouncing ball orbit between the components is determined by 
its resonances (poles of its scattering matrix) among other such exterior 
domains [ZIO]. The proof is essentially the same as in the interior case, 
once some known facts on resonance poles (explained to the author by 
M. Zworski) are added. 

4.2.1 Domains and metrics with the same wave invariants 

A Penrose mushroom type example due to Michael Lifshits (see [Me] , [Rau] 
for pictures and background) shows that wave invariants arc not sufficient 
to discriminate between all pairs of smooth billiard tables. Indeed, Lifshitz 
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constructs (many) pairs of smooth domains (Qi,Jl2) which have the same 
length spectra and the same wave invariants at corresponding pairs of closed 
billiard orbits 7j of Qj (j = 1,2). It follows that the Poincare maps "y^. have 
the same Birkhoff normal forms as well. 

The idea is to exploit the complete integrability of the biUiard flow on 
an ellipse, i.e the fact that it is foliated by caustics. Caustics are curves in 
the domain with the property that any billiard trajectory which starts off 
tangent to the caustic will remain tangent to it, as with the confocal ellipses 
and hyperbolae of an ellipse. The billiard trajectories fall into two families 
separated by the bouncing ball orbit between the foci. One family consists 
of trajectories which remain tangent to confocal ellipses (which degenerate 
to the segment between the foci). The second family consists of trajectories 
which are tangent to confocal hyperbolae that pass through the segment 
between the foci. Now divide the major axis into two parts, the segment 
between the foci and the other two segments. Remove the segment between 
the foci and replace it with any simple curve with the foci as endpoints, 
e.g. a 'tongue' below the segment. Next, remove the outer segments of the 
axis between the foci and replace them with any 'bumps' below the segment. 
The trajectories which start in the outer bumps never intersect the segment 
between the foci and therefore never go into the tongue. Similarly trajectories 
which come into the elliptical part from the tongue pass through the segment 
between the foci and never go into the outer bumps. 

It follows that the closed billiard orbits fall into two families: those which 
never intersect the segment between the foci and therefore bounce back and 
forth between the bumps; or those which do intersect this segment and never 
go into the outer bumps. 

To obtain non-isometric domains with the same wave invariants it suffices 
to reverse the relative orientations of the two families, either by reflecting the 
segment between the foci at the center (i.e. mirror reversing the tongue) or 
equivalent ly by reversing the outer bumps. Since the reversal involution 
is an isometry on each 'half of the domain, it does not change the wave 
invariants for each half. But there is no 'interaction' between the halves (i.e. 
no closed geodesic intersecting both halves), so the involution preserves all 
wave invariants. We refer to [Me, Rau] for pictures. 

The question may occur whether such a domain and its reversal are 
isospectral or not. This seems dubious, but we are not aware of a proof of it. 
As mentioned above, equality of wave invariants implies equality of Birkhoff 
normal form invariants at corresponding pairs. Are the billiard maps of each 
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pair symplectically conjugate? Again dubious, but again we don't know a 
proof. 

4.3 Isospectral deformations and spectral rigidity 

One of the first uses of wave invariants and dynamics was the proof by 
Guillemin-Kazhdan [GK] that negatively curved surfaces are spectrally rigid. 
In the negatively curved case, the Maslov indices arc always zero and no 
cancellation takes place in the wave trace formula as one sums over closed 
geodesies of the same length. Hence, Lsp{M,g) is a spectral invariant of 
negatively curved manifolds. An isospectral deformation therefore preserves 
the length spectrum. 

If 7 is an isolated, non-degenerate closed geodesic of then for any 
deformation gt of g, 7 deforms smoothly as a closed geodesic 74 of gt and one 
may define its variation 

= j^\t=oL^,. (4.9) 

It is not hard to compute that 




(4.10) 



where g is viewed as a quadratic function on TM and 7 is viewed as the 
curve (7(s),7'(s)) in TM. 

It follows that whenever the closed geodesies are non-degenerate and of 
multiplicity one in the length spectrum, we have 

gds = 0, V7. (4.11) 

Thus, the Livsic cohomology problem enters. Since Lsp{M, g) is generically 
simple, this equation holds for generic isospectral deformations. We note 
that only the principal level of the wave invariants was used here. 

In the case of negatively curved surfaces, the geodesic flow is Anosov and 
it is known that the cohomology is trivial, i.e. that g = S(/) for some smooth 
/ (the regularity problem is a separate issue, but for the sake of brevity we 
will not consider it). The next step is to study harmonic analysis on the unit 
sphere bundle S*M to determine if there actually can exist / when g is a. 
quadratic form. For surfaces of negative curvature, it was proved in [GK] that 
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there cannot exist a smooth solution, and hence there exist no isospectral 
deformations of negatively curved surfaces. After a series of partial results by 
Guillemin-Kazhdan, Min-Oo and others, the result was extended to higher 
dimensions by Croke-Sharafutdinov [CS]. In [SU], Sharafutdinov- Uhlmann 
prove infinitesimal spectral rigidity for the more general class of closed two- 
dimensional manifold without focal points whose geodesic flow is of Anosov 
type. 

Isospectral deformations also lead to the inverse marked length spectral 
problem and hence (in some cases) to the conjugacy rigidity problem for 
the geodesic flows. An isospectral deformation preserves lengths of closed 
geodesies in the generic case where the multiplicities all equal one (or more 
generally where no length is cancelled in the wave trace formula), and there- 
fore marks the length spectrum, i.e. gives a correspondence between closed 
geodesies and lengths, by the length spectrum of the initial metric. More 
precisely, as one deforms the metric or domain, each one-parameter family 
of closed geodesies % stays in a fixed free homotopy class of the funda- 
mental group of M. Therefore, an isospectral deformation gives rise to a 
one-parameter family of metrics or domains with the same marked length 
spectrum. As discussed above, the marked length spectrum determines the 
metric for surfaces of negative curvature (Croke, Otal), which gave a different 
proof of the result of [GK]. 

There do not seem to exist many studies of the Livsic cohomology equa- 
tion for non-hyperbolic systems. It was proved by Kuwabara using heat 
invariants that fiat tori are spectrally rigid. It might be interesting to review 
the result in terms of the Livsic equation, and possibly generalizing them 
to a broader class of completely integrable systems. As will be discussed in 
§8.1, closed geodesies of integrable systems (such as surfaces of revolution or 
Liouville tori) come in families on invariant tori T, and the analogue of (4.11) 
would say that 



for each invariant torus of the geodesic flow. This is an analogue of the 
Livsic cohomology equation when geodesies come in non-degenerate criti- 
cal manifolds. In the case of integrable systems, it could be studied using 
Fourier analysis relative to the foliation by invariant tori of a completely in- 
tegrable system. However, it is weaker than the Livsic equation. At least 
for some Riemannian surfaces, isospectral deformations of metrics with in- 
tegrable geodesic flows must be through metrics whose geodesic flows are 




(4.12) 
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C° integrable (see §8.1), and the study of the Livsic type equation (4.12) 
might combine to form a more powerful tool. The simplest case to consider 
is whether simple surfaces of revolution on S'^ are spectrally rigid (see §8.1 
for background). 

To give an extreme example of an open isospectral deformation prob- 
lem, consider the case of ZoU manifolds. It is simple to see that isospectral 
deformations of Zoll manifolds must be ZoU [Z8], but it is an open problem 
whether any non-trivial isospectral deformations exist even for M = S'^. One 
of the main problems is that all principal symbol level spectral invariants of 
Zoll manifolds are the same, so one has to dig further into the wave trace 
expansions to find obstructions to isospectral deformability. 

4.3.1 Isospectral deformations of domains 

A few words on the analogues for Euclidean domains. To the author's knowl- 
edge, there are no spectral rigidity results even for Euclidean plane domains 
except in the case of a disc. It is not even known whether ellipses are spec- 
trally rigid. The length spectrum of a bounded domain is a spectral invariant 
for both Dirichlet and Neumann boundary conditions, so an isospectral de- 
formation must preserve the length spectrum and even the marked length 
spectrum. We will discuss some results of K.-F. Siburg on the marked length 
spectrum problem in §8.4. 

Since the length spectrum is generically an isospectral invariant, it is 
natural to consider the analogue of the Guillemin-Kazhdan rigidity result 
for domains with hyperbolic billiards. The billiard flow is hyperbolic, so 
presumably the Livsic equation can be solved (although the author is not 
aware of a specific reference). It is not clear how the harmonic analysis 
would change from the case of negatively curved surfaces. 

There is another approach which reduces the wave group to the boundary. 
Suppose that fit is a smooth one parameter family of smooth compact plane 
domains such that the Dirichlet (resp. Neumann) spectrum Spec{A^'^^) 
(resp. Spec{A^^^^) ) is constant in t. Since the area and length of the 
boundary are spectral invariants, both must be fixed under the deformation. 
With no loss of generality, we assume that the variation is generated by a 
normal vector field pu, where i/ is the outward unit normal to Q(0) and where 
pe C~(9fi(0)). 

The variations of the eigenvalues are given by Hadamard's variational 
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formulae [0]: 



Ian(o) P \du(i)j{0)\dn\^dA, Dirichlet 
A,(0) = <( (4.13) 
IdnP l</'j(0)|9f^(o)P dA, Neumann. 

Hence, the infinitesimal deformation condition is that the right hand sides 
are zero for all j. To normalize the problem, we assume with no loss of 
generality that the deformation is volume preserving, which implies that 

pdA = 0. (4.14) 

an 

Any such p defines a volume preserving deformation of Q. 

Thus, the infinitesimal deformation is orthogonal to all boundary traces 
of eigenf unctions: 

{/an(o) P \du4>j{0)\dn\'^dA ^ 0, Dirichlet 
(4.15) 
lanP l^j(0)l9f2(o)P dA = 0, Neumann. 

We may rewrite these conditions in terms of the boundary values of the wave 
kernel: 

{du^d,yyU{t,x,x)\a;^dn, Dirichlet 
(4.16) 
C/(i, X, x)|a;gaQ, Neumann 

as saying that 

E^t, X, x)p(x)dA(x) = 0, yt. (4.17) 



'an 

Using the calculations in [0], one can obtain expansions for (4.17) at 
t — similar to a heat kernel expansion in terms of integrals of p against 
polynomials in the (extrinsic) curvature invariants of the boundary. The 
singularities at i 7^ in turn give integrals over closed orbits of the billiard 
map. In the case of an ellipse, one can directly study the boundary traces of 
the eigenfunctions, and it appears that Jp pds = for every caustic T for the 
billiard map. It would be interesting to explore this further. 



40 



5 Formulae for wave invariants 



When studying the inverse spectral problem for individual pairs of metrics 
rather than isospectral deformations, it is difficult in general to bring in 
the global dynamics of the geodesic flow. Nothing better is known than to 
concentrate on a single closed geodesic. So to make progress on the inverse 
spectral problem, it is crucial to have useful formulae for the wave invariants 
associated to a closed geodesic 7 and its iterates Y- The purpose of 
this section is to survey the known methods of calculation and the formulae 
which they bring. 

There are several potential approaches: 

• (i) Construct a microlocal parametrix for e'*^ at 7 and apply a sta- 
tionary phase method to calculate the wave invariants (cf. [Z9]). 

• (ii) Construct a Birkhoff normal form for e**^ at 7 and relate normal 
form invariants to wave invariants (cf. [G, Zl, Z3, Z4, Z9]. 

• (iii) Construct a Birkhoff normal form for the monodromy operator M 
on the microlocal solution space ker^(j(A — A^) at a fixed initial point 
mo of 7 and relate it to the wave invariants (cf. [SjZ, ISZ]). 

• (iv) For bounded domains, apply the Balian-Bloch or Calderon projec- 
tor methods (cf. [Z5]§7). 

Wc now display the formulae and give a brief discussion of the methods. 
In the following sections, we describe the methods in much more detail. 

5.1 The parametrix method 

In boundaryless manifolds one can construct a Hadamard parametrix for 
cos t\/A{x, x) for small times. This operator solves the initial value problem 




(5.1) 



and has the form of the real part of the oscillatory integral 




mod C" 



100 



(5.2) 
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where the Hadamard-Riesz coefficients Uj are determined inductively by the 
transport equations 

e' TT I dUo _ n 
20^0 + — - U 

(5.3) 

The solutions are given by: 

Uo{x,y) = e-'/\x,y) 

(5.4) 

C/,+i(x, y) = e-'/\x, y) s^e{x, x^f'^A^U^ix, x,)ds 

where Xg is the geodesic from a; to y parametrized proportionately to arc- 
length, where Q{x, y) denotes the volume density in normal coordinates cen- 
tered at X, and where A2 operates in the second variable. 

The simplest case is where the metric is without conjugate points, in 
which case the parametrix is global in time on the universal cover. It may 
then be projected to M in the usual way by summing over the deck trans- 
formation group. One then takes the trace and computes the coefficients by 
a stationary phase method. 

This was first done by Donnelly [D] for negatively curved surfaces and 
then for all manifolds without conjugate points by the author in [Z9]. To 
get an impression of the complexity of the result, here is the formula for the 
subprincipal wave invariant (the coefficient of the logarithmic singularity) : 

+ Hess(^);i(e;^((0,0),((j,0)d(j 

+ C"n,o,o,l t t («<70{HeSs(/,);l}2(^,e; V[,]) i=,ds 

+ Jo"^ e|(ai)(A2e-i)((ai, 0), (a2, 0)) • da.da^] . 

Here, Q^{x) = 9(a;, 7a;), is a similar quantity using geodesic polar coordi- 
nates, r = r{x, y) denote the distance function of the universal cover, 7 G F, 



42 



the fundamental group, /^(x) = r{x, denotes the displacement function. 
Sec [Z9] for details. The higher wave invariants are of course exponentially 
more complicated. 

In the case of bounded domains, the parametrix approach is even messier. 
A microlocal parametrix for cost-\/A near a transversal reflecting ray was 
constructed by J. Chazarain [Ch] (see also [GM, PS] for more details ). In 
principle it could be used in the same way as a microlocal parametrix in 
the boundaryless case, but in the boundary case there are almost always 
conjugate points and additionally the transport equations and phase are 
more complicated. Applying the method of stationary phase to all orders 
to such a parametrix descends into a jungle of formulae. To our knowledge, 
no concrete inverse spectral results have been proved using the parametrix 
method. 

5.2 Wave invariants and quantum Birkhoff normal 
forms 

Evidently, some guiding principle is needed to civilize the wilderness of for- 
mulae. Two such principles have emerged, at least in the case of bounded 
domains. The first is the method of quantum Birkhoff normal forms. As 
for classical Birkhoff normal forms, one conjugates the Laplacian to a model 
normal form on a model space. The wave invariants are then expressed in 
terms of the coefficients of the normal form. The method is described in 
detail in the next section. For the moment, we are interested in the formula 
it gives for the wave invariants. 

To state the result we will need some notation. The wave invariants at a 
closed geodesic 7 are invariants of the germ of the metric at 7. For simplicity, 
we assume that the geodesic is isolated and non- degenerate (det (J — P-y) ^ 0). 
We assume dimM = n + 1 and introduce Fermi normal coordinates {s^y) 
along 7. We denote the corresponding metric coefficients by Qij, and refer to 
the coordinate vector fields ^ (j = 1, . . . , n) and their real linear com- 
binations as Fermi normal vector fields along 7. We denote the Riemannian 
connection, rcsp. curvature tensor, by V resp. R and refer to contractions of 
tensor products of the V™-R's with the Fermi normal vector fields as Fermi 
curvature polynomials. Such polynomials will be called invariant if they are 
invariant under the action of 0(n) in the normal spaces. 

We consider the Jacobi equation Y" -\- R{Y, ^')^' — along 7. The space 
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of complex Jacobi fields along 7 is denoted The linear Poincare map is the 
monodromy map Y(t) Y(t + Ly) on this space. We refer to its eigenvectors 
Yj, Yj as Jacobi eigenvectors. We denote by yjk their components relative to 
the Fermi normal vector fields. 

We then define Fermi-Jacobi polynomials to be invariant contractions of 
tensor products of V"*-R's against ^ and against the Jacobi eigenvectors, 
with coefficients given by invariant polynomials in the components yj^. We 
will also use this term repeated integrals over 7 of such polynomials. Finally, 
FJ polynomials whose coefficients are given by polynomials in the Floquet 
invariants 

are Fermi-Jacobi-Floquet polynomials. Here, we assume for simplicity that 
the closed geodesic is elliptic and that the eigenvalues of its Poincare map 
are {e*"^}. (The Floquet invariants in the general non-degenerate case are 
analogous; see [Z4] for details). 

We define the 'weight' of the data going into a Fermi-Jacobi-Floquet poly- 
nomial in terms of its scaling behavior under g — > e^g. Thus, the variables 
gij, D^ygij (with m := (mi, . . . , nin+i)), L := L^, aj, ytj, ytj have the follow- 
ing weights: wgt{D'^ygij) = -|m|, wgt(L) = 1, wgt(aj) = 0, wgt{yij) = 
|, wgt(yjj) = — |. A polynomial in this data is homogeneous of weight s if all 
its monomials have weight s under this scaling. 

The general result is [Z4]: 

Theorem 5.1 Let 'y be a strongly non-degenerate closed geodesic. Then 

a^k = IrA^'^i 9)ds where: 

(i) I^.k{s;g) is a hom,ogeneous Fermi-Jacobi-Floquet polynomial of weight 
—k — 1 in the data {y^j, ij^j, D"^yg} with m — (mi, . . . , m„+i) satisfying \m\ < 
2k + A ; 

(a) The degree of L/-k in the Jacobi field components is at most Qk + 6; 
(Hi) At most 2k + 1 indefinite integrations over 7 occur in I^-k,' 
(iv) The degree of I^-k in the Floquet invariants (3j is at most k + 2. 

The formula is simplest in dimension 2, where there is only one Floquet 
invariant (3. We use the notation r for the scalar curvature, Ti, for its unit 
normal derivative along 7, t^i^ for the Hessian Hess(r)(z/, z/). We denote 
by Y the unique normalized Jacobi eigenvector along 7 and by Y its time- 
derivative. 
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The subprincipal wave invariant a^o is then given by: 



j^o — 



L* 



(5.5) 



where: 

(a) a-y _i is the principal wave invariant (4.8); 

(b) L# is the primitive length of 7; a is its Morse index; Pj is its Poincare 

map; 

(c) BjQ.j has the form: 



B. 



■70;j 



+ 



1 r^* 

/ [a |F|^+6i T\Y-Y\'^+b2 rRe {YYf+c T^\Y\'^+d r^^|r|^+e SjoT]ds 
L* Jo 



L* 



E 



Cl-r. 



sin((n 



mja] 



0<m,n<3;m+n=3 



:,i{m—n)a^ |2 



r^(s)y"^ • y'^(s)ds 



+ 



E 



0<m,n<3;m+n=3 




ds 



for various universal coefficients. 

The formula is very complicated, and it is only the second term in the 
expansion! That is why one hopes to directly use the information that these 
coefficients determine the Birkhoff normal form. 

In the case of surfaces of revolution, it is possible to simplify the formulae 
to the point where one can determine the metric along the invariant geodesies. 
The point is to recover the Taylor expansion of the curvature and its normal 
derivatives along invariant geodesies. In [Z2], a variant of this method was 
used to show that 'simple' surfaces analytic of revolution (those with just 
one invariant geodesic) are determined by their wave invariants, hence are 
spectrally determined among such surfaces. The variant was to use the exis- 
tence of a global Birkhoff normal form for the metric and Laplacian, which 
greatly simplified the calculations. 

The wave invariants can also be calculated by the method of Birkhoff 
normal forms in the boundary case, but the method becomes appreciably 
more difficult. We therefore use a different approach which is inspired by the 
Balian-Bloch approach to the Poisson formula. 



45 



5.3 Balian-Bloch approach in the boundary case 



We now consider wave invariants associated to periodic reflecting rays of a 
bounded plane domain Q. For simplicity and because they are often useful 
in applications, we consider wave invariants at a a 2-link periodic reflecting 
ray 7 of length rL^ (a 'bouncing ball orbit'). We refer to the §3.3 for the 
basic definitions. 

We introduce some further notation: we write y = fa{j){xj) where a{2j + 
1) = — , — +. The length functional in Cartesian coordinates for a 

given assignment cr of signs is given by 



2r-l 



. . . , X2r) {Xj+l - XjY + {fa{j+l){Xj+l) - fa{j){Xj)f. (5.6) 

i=i 

We will need formulae for the entries of its Hessian in Cartesian coordi- 
nates at the critical point corresponding to the rth repetition of the bouncing 
ball orbit with a given orientation. Thus, Xj — for all j. We will assume, 
with no essential loss of generality, that q{xodd) — -A, q{Xev) — B. We put: 
a = l-L/^(0),6 = l-L/!^(0). 

The Hessian is given in either angular or Cartesian coordinates by: 

a 1 ... 1 



H2r — — \ 



(5.7) 



where a 



2(1-^)' 



1 b 1 ... 

1 a 1 

16 1... 

1 ... b 



2(1 — ■^). A well known formula relates the 



determinant of this Hessian to that of the Poincare map: 
det(/ - P^r) = L^'- det H2r; (7 2 - link.) 



(5.8) 



We denote the matrix elements of the inverse Hessian if2r^ ^^e bouncing 
ball orbit by h^^ . To be precise, the bouncing ball orbit has two possible 
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orientations, one (which we denote by + which starts at the top graph y — 
and proceeds to the bottom y = f- and the other — which in the reverse 
order. We then have two (closely related) length functions and Hessians. We 
denote the matrix elements of their inverses by h^J.. One has = hF^r^''^~^. 
We now state the formulae for the wave invariants. 

Theorem 5.2 [Z5] Let j be a 2-lmk reflecting ray of length L^, and let Y) 
resp. 7"'' be the rth iterate of 'j, resp. •y^'^. Then there exist polynomials 
P2,r,jiCi, ■ ■ ■ , C2j', ^1, ■ ■ ■ , V2j) which are homogeneous of degree —j under the 

dilation f A/, which are invariant under the substitutions <^=^ rij and 
under f{x) f{—x), and which have degree j + 1 in the Floquet data e^"^, 
such that 

«r J = P2,rAf-\Q): /- ^0), • • • , f^y^'\oy, /f (0), /f (0), • • • , /f +'^(0)). 

The leading order term in derivatives of f± has the form 

r{{hfy EjLi;,.o hi' + {hfy-' E:,.o(^?)n/f "'^(o)/f (0) 

+rmy E:,.i hi' + (hiy-' El.^{hlr}f^^^''\0)f^^\0) 
-r{{hfy-'hl' Y.Un^.h'' + {h'^y-'T.U-.^Ahmft'-'\^)f'^\^) 

-rmy-'h^:Y.Un^oh'i + {h'^r'Y.U^^^^ 

+R2r{f'-'f+iO),f^-'f40)), 

where the remainder R2r{j^-'~^f+{0), j^-'~'^f-{0)) is a polynomial in the des- 
ignated jet of f±. 

In the case where the domain has a symmetry interchanging the top and 
bottom, so that f± are mirror images of the graph oi y = f{x). The length 
functions for Y ^^'^ 7 ' (and hence their Hessians) are equal, so we may 
drop the subscripts ±. The formula simphfies as follows: 

Corollary 5.3 Suppose that y (as above) is invariant under an isometric 
involution a. Then, modulo the error term R2rij^-'~^fiO)), we have: 

a,r,_, = r{2{h'y f^^^KQ) + {2{h''y.,^^, + {h''y-'Y.U{h''f}f''\^)f''-'\m}- 

Here, we sum over repeated indices. 
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The inverse spectral problem is then reduced to analyzing Hessian coef- 
ficients. It is not evident to the author how this reduction of the inverse 
spectral problem would be visible using Birkhoff normal forms. 

5.3.1 Melrose-Marvizi invariants 

In [MM], Melrose-Marvizi introduce further spectral invariants of a con- 
vex domain which could be interpreted as quantum normal form invariants, 
namely the normal form of A around the closed geodesic dQ. The calcula- 
tion in [MM] involves Melrose's normal form for glancing hypersurfaces and 
may be viewed as a normal form construction, but one might also follow 
Lazutkin's construction of whispering gallery quasi-modes to put A into nor- 
mal form as a function of a so-called Airy operator. This would presumably 
give a new way to calculate the invariants. J. Toth and the author have 
partial results in this direction (unpublished), and from discussions with G. 
Popov it appears that they are known to others. 

The first two Melrose-Marvizi integrals have the form [MM]: 



6 Calculation of wave invariants I: Birkhoff 
normal forms 

We now outline the construction of the Birkhoff normal form of A in both the 
boundaryless and boundary cases, and indicate how it is used to calculate the 
wave invariants. Our goals are to explain how the calculations in Theorem 
5.1 are done. 

Let 7 be a non-degenerate closed geodesic on an n-dimensional Rieman- 
nian manifold, and at first let us assume it to be elliptic. During the 70's, 
various authors (Babic, Lazutkin, Ralston, Guillemin-Weinstein) constructed 
a series of quasi-modes and quasi-eigenvalues associated to 7. Roughly speak- 
ing, the results showed that for each transversal quantum number q £ Z"'~^, 
there was an approximate eigenvalue of the form 




2 



+ 



(6.1) 
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where 




where the coefficients are polynomials of specified degrees and parities. We 
refer to [BB] for a clear exposition and for details. 

A natural question is whether the wave invariants a^j of (4.6) can be de- 
termined from the quasi-eigenvalues (6.1), i.e. from the polynomials Pj{q), 
which could be computed explicitly from the quasimode constructions de- 
scribed in [BB]. This question was difficult to answer or even to formulate 
precisely until the article of V. Guillemin [G] on quantum Birkhoff normal 
forms for the Laplacian around elliptic closed orbits. This article does not 
refer to quasi-modes or quasi-eigenvalues but it effectively proves that the 
wave invariants may indeed be expressed in terms of the polynomials Pj{q) 
and conversely that these polynomials are spectral invariants. A somewhat 
different proof of this result, which constructed quantum Birkhoff normal 
forms by adapting the algorithm in [BB] for constructing quasi-modes, was 
given by the author for elliptic orbits in [Z2], and for general orbits in [Z3]. 

The results of [G] are stated in the terms of quantum Birkhoff normal 
forms. A number of expositions now exist which explain this notion in detail 
(see for instance [Z9] in addition to the original articles), so we will only 
briefly review the notion. To put A into normal form, is ffist to conjugate 
it into a distinguished maximal abelian algebra A of pseudodifferential op- 
erators on a model space, the cylinder S"], x R", where S}^ is the circle of 
length L. The algebra is generated by the tangential operator Dg := 
on S}^ together with the transverse action operators. The nature of these 
action operators depends on 7. When 7 is elliptic, the action operators are 
harmonic oscillators 



When 7 is non-degenerate, they involve some mixture of these operators (and 
also complex hyperbolic actions) according to the spectral decomposition of 
Pj. For notational simplicity we restrict to the elliptic case and put 




while in the real hyperbolic case they have the form 




k=l 
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where e^'"'= are the eigenvalues of the Poincare map P^. 

To put A into normal form is to conjugate it to the model space and 
algebra as a function of Dg and the action operators. The conjugation is only 
defined in a neighborhood of 7 in T*M — 0, i.e. one constructs a microlocally 
elliptic Fourier Integral operator W from the conic neighborhood of R+7 in 
T*N^ - to a conic neighborhood of T*Sl in T*{Sl x RJ") such that: 

.-1 _ m , ■■■Jn) , P2(A, ■■■Jn), , Pk+l{h, ■■■Jn) 



Ub + (LD)2 (L2))^+i 

(6-2) 

where the numerators In),Pj{h, ■■■,in) are polynomials of degree j+1 

in the variables /i, /„, where denotes a microlocal inverse to W. Here, 
2) = Ds + j^Ha. The kth remainder term lies in the space <S)'^t'^02{k+2-j)'^^~-' , 
where denotes the pseudo-differential operators on the model space of 
order s and where 02(k+2-j) denotes the operators whose symbols vanish to 
the order 2{k + 2 — j) along 7. We observe that (6.2) is an operator version 
of (6.1). 

The inverse result of [G], see also [Z3, Z4]] is: 

Theorem Let be a non-degenerate closed geodesic. Then the quantum 
Birkhoff normal form around j is a spectral invariant; in particular the clas- 
sical Birhoff normal form is a spectral invariant. 

In other words, one can determine the polynomials Pj{Ii, .... In) from the 
wave trace invariants of A at 7. When 7 is a non-elliptic (e.g. hyperbolic) 
geodesies, the action operators have continuous spectra and there are no 
approximate eigenvalue expansions as in (6.1). Thus the quantum Birkhoff 
normal form approach is conceptually clearer. 

We give a brief review of how wave invariants and normal form invariants 
for the Laplacian are calculated by the algorithm in [Z2, Z3, Z4, Z5] in the 
case of an elhptic closed geodesic. The algorithm is similar in the general 
non-degenerate case. We then consider the algorithm for calculating the 
normal form of the monodromy operator in [ISZ]. 

6.1 Manifolds without boundary 

Our algorithm for conjugating A to normal form along a closed geodesic 7 is 
inspired by the constructions due to Lazutkin [L] and Babich-Buldyrev [BB] 
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for constructing quasimodes associated to elliptic closed orbits. The same 
conjugation method works in the general non-degenerate case. 

Since we have recently written an exposition of the method [Z9], we only 
include only a few formal aspects of the calculations, in the hope they provide 
sufficient explanation of how the wave invariants are calculated. In particular, 
we wish to emphasize: 

• The semiclassical normal form, and how it arises from a a semi-classical 
scaling of the Laplacian along 7. Our method differs from the others 
[G, ISZ] by working entirely on the quantum level and inductively on 
the Taylor expansion of the metric around 7. The homogeneous normal 
form (6.2) is obtained from the semiclassical normal form by (roughly 
speaking) replacing the large parameter by \Ds\ along 7. 

• The homological equations and the obstructions to their solvability 
along a closed orbit. The obstructions determine the normal form. 

6.2 Quasi- mode heuristics 

We wish to conjugate the Laplacian to a normal form on a model space, 
namely the normal bundle of 7, or equivalently, the cyhnder 5"^ x M", 
where as above Sj^ = M/LZ. On the phase space level, the model is T*[S\ x 
M") or more precisely the cone \ri\ < to, \y\ < e in the natural symplectic 
coordinates (s, cr, y, r\) corresponding to Fermi normal coordinates along 7. 

In the construction of quasi-modes associated to 7, one uses the WKB 
ansatz 

^kq{s, ^/¥^y) = e'''>"^'Ukg{s, ^/¥^y, r^^), 
where rkq — ^{k + {q + ^)a) and where 



00 




j=0 



One then solves asymptotically the eigenvalue 

Aye^'^'^^'Ukgis, ^qy, r^^) ~ \kqe'''*"''Ukq{s, ^qy, r^^), 

satisfying the periodicity condition of being well defined on the cylinder. The 
intertwining operator W^^ to normal form may be thought of as the operator 
taking the model eigenfunctions to the quasimodes, 
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We now change our point of view and concentrate on the construction of the 
intertwining operator. 



6.3 Semiclassically scaled Laplacian 

In view of the form of the quasi-eigenvalue problem, 

i 1 i 1 

A^ChL'^U {s, h~2u, h) = \{h)e'hL''U (s, /i" 2^, h), 

it is natural to rescale the Laplacian before conjugating it to normal form. 
In fact, this rescaling is all we retain from the quasi- mode construction. We 
therefore introduce the unitary operators Th and Mh on !Kt or equivalently 
on the 1/2-density version L|,(M-^ x M"',r2i/2) given by 

(2.5a) Th{f{s,u)\ds\^'''\du\^''') -.^ h'''/'' f{s,h--^u)\ds\^'''\du\^''' 



(2.56). Mh{f{s,u)\ds\^'^\du\^'^) ■.= e^'f{s,y)\ds\^'^\du\^l^ 

Definition The rescaling of an operator = a{s, Dg, u, Dy) of the adapted 
model is the operator 

Ah := T*M*AnMh. 



We now rescale the Laplacian in Fermi normal coordinates. It is conve- 
nient to first conjugate to the unitarily equivalent 1/2-density Laplacian 

which can be written in the form: 

n 

(2.9) - Ai/2 = J- V^a./Va, J- V2 + J2 J-''^du,g'^Jdy. J" 

ii=l 

n n 
ij=l 1=1 

Here, ■= and J = J(s,u) = ^ is the volume density in these coordi- 
nates. 
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We then have: 

(2.10) -M*AMh = -{hLY^g"" + 2i{hL)-^g°''ds + i(/iL)-^r° + A 

Conjugation with then gives 
(2.11) 

n n 
ij=l i=l 

the subscript [h] indicating to dilate the coefficients of the operator in the 
form, fh{s,u) := f{s,h^u). 

Expanding the coefficients in Taylor series at h — 0, we obtain the asymp- 
totic expansion 

oo 

(2.12) A, ~ 5^ 

m=0 

where £.2 = L~^, £3/2 — and where 

j=l ij=l 

6.4 Conjugation of scaled Laplacian to semi-classical 
normal form 

We now conjugate the semi-classically scaled Laplacian to normal form. The 
conjugating operators are s-dependent operators acting on the transverse 
space to 7 (they are sometimes called tangential semi-classical Fourier inte- 
gral operators). 

The first step is to conjugate the principal term into quadratic normal 
form on the transverse space. Since it is quadratic, there exists a metaplectic 
conjugation ii{s) depending on s on the transverse space which conjugates it 
to the operator D of (6.2). That is, for each s e [0, L] there exists an operator 
n{s) in the metaplectic representation of the metaplectic group ML{n, M) — > 
Sp{n, M) which acts on the transverse of the model space. For background 
on metaplectic operators and details on the conjugation we refer to [Z3, Z4]. 
This conjugates the other operators Jij to new operators D^. 
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Once the principal term is in normal form, we continue by conjugating the 
lower order terms to functions of the actions by using perturbation theory, 
i.e. formal series in h. Thus we wish to put the formal series Ah into a 
semi- classical normal form by an infinite sequence of conjugations, each one 
putting one new order in h into normal form. We carry out the procedure to 
two orders to illustrate the main points. 

The first step is to construct Q i (s, D^) such that 

^-..ig, ^^^..ig, |„ = [-h-^ + 2h-'V + 2)1 + . . . ]|„ 

where 

|o denotes the restriction of the operator to functions of y, (6.3) 

and where the dots . . . indicate higher powers in h. Introduction of the |o 
operator is motivated by the construction of quasi-modes, since the normal 
form is only being applied to the amplitude (a function oi y). A more con- 
ceptual explanation (suggested by the work of lantchenko-Sjostrand-Zworski 
[ISZ]) is that the correct Hilbert space on which to define the operators is 
the microlocal solution space of A — along 7. These are essentially the 
space of quasi-modes along the open arc (0, L) of 7. Thus, |o restricts the 
operator to the microlocal solution space. 

Expanding the exponential, we find that the operator Qi then must sat- 
isfy the homological equation 

{[L-^2),Qi] + Di}|„ = 0. 

One may solve this equation explicitly by further conjugating D to Dg. When 
T) = Dg, the homological equation becomes 

{[L-'Ds,Qi]+'Di}l^O, 

that is. 

Here, dgA is the Weyl operator whose complete symbol is the s-derivative of 
that of A. To solve this equation we rewrite it in terms of complete Weyl 
symbols. We will use the notation A{s, x, ^) for the complete Weyl symbol 
of the operator A{s, x, Dj.). We then arrive at the homological equation 

L-^dsQ^{s,x,i) = -iDi|o(s,x,0 
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We solve with the Weyl symbol 

Qi{s,x,^) ^ Qi{0,x,^) + L / -i'Di\o{u,x,^)du 
^ ^ Jo ^ 

where Qi{0,x,^) is determined by the consistency condition 

Qi{L,x,i) - Qi{0,x,i) = L / -i'Di\o{u,x,i)du.) 
^ ^ Jo 

To solve the equation, we invoke the fact (which is not obvious) that 
T)i\o{u, X, ^) is a polynomial of degree 3 in {x, We also switch to complex 
coordinates Zj = Xj + i^j and Zj = Xj — i^j in which the action of ra{L) is 
diagonal. The homological equation becomes 

Qi{0,e'''z,e-'''z)-Qi{0,z,z)=L [ -I'D i\o{u, z, z)du. 



We put: 



2 ' 

|m|+|n|<3 



and 

|m| + |n|<3 

then the homological equation becomes 



|m| + |n|<3 |m| + |n|<3 



Since there are no terms with m = n in this (odd-index) equation, and since 
the a/s are independent of tt over Z, there is no obstruction to the solution. 
Thus we can simply eliminate the term Di. 

6.4.1 The normal form coefficients 

We now consider the second step, where the first obstruction occurs. We thus 
seek a pseudodifferential operator Qi{s,x,Dx) and a quadratic polynomial 
/o(/i, /„) in the action operators so that 

2 
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with 

As usual, the dots signify terms of higher order in h. The homological equa- 
tion is then 

{p„gi] + D|}|, = /„(/!,...,/„), 

or equivalently 

dsQl\o = {-Di +fo{h,...,In)}\o. 

We rewrite the equation in terms of the complete Weyl symbols, to obtain 

1 

L~^dsQl{s, Z, Z) = \o{s, Z, Z) - /o(|^l|^ ■ ■ ■ , l^nH} 

or equivalently 

Qi{s,z,z) = Qi{0,z,z) -iL / ['Di\o{u,z,z) - fo{\zi\'^, . . . ,\zn\^)]du, 

Jo 

or again, 

1 

gi(0, e'^'z, e-*"^)-gi(0, z, z) = -iL{ / 2)| z, z)du-Lfo{\z,\^ . . . , \znf)}. 

Jo 

1 

By construction. Do \o{u, z, z) is a polynomial of degree 4, so if we put 

Ql{s,Z,z)= qi;mn{s)z'^z'', . . . , = ^ C^fcl^l^^^ 

|m|+|n|<4 |fc|<2 

and 

1 _ 1 _ _i 1 I 

'DS\o{s,Z,z)du -.^ dS-mn{s)z'^z'', dS-mn '-^ dS-mn{s)ds, 

|m| + |n|<4 

we can solve for the off-diagonal coefficients. 



We also must set the diagonal coefficients equal to zero, and this determines 
the Cok coefficients: 

Cok = '^l;fcfe- 

These coefficients are the quantum normal form coefficients. Since we used 
only algebraic operations on the rescaled Laplacian, it is clear that the coef- 
ficients can be calculated in terms of the data described in Theorem 5.1. We 
refer to [Z4, Z9] for further details. 
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6.5 Bounded domains 

There are two serious differences in tfie normal form for the Laplacian around 
a bouncing ball orbit of a bounded domain. First, we need to straighten out 
the domain to define a model space. Secondly, we need the conjugation to 
normal form to incorporate the boundary conditions. 

6.5.1 The model domain Qo 

The configuration space of the model is the infinite strip fig. We denote the 
coordinate on [0, L] by s and that on M by |/ , with dual cotangent coordinates 
a,r} on T*[0,L] x T*'R. As in the boundaryless case, we would like to view 
the model space as the normal bundle of the orbit. In the case of a bouncing 
ball orbit in the boundary case, the normal bundle and exponential map 
are ill-defined at the reflection points but Lazutkin has constructed a nice 
replacement for them. Namely, he constructs a formal power series map from 
fig to Qq, 

^-.n.^Ue, ^s,y)^{s,y) 

of the form 

y = E^l^pp(s)y^ 

with real valued analytic coefficients in s extending analytically to a neigh- 
borhood of [0, L] in C and satisfying the boundary conditions: 

{s = 0}U{s = L} = $(aa). 

Lemma 6.1 Suppose that AB is a bouncing hall orbit. Then there exists a 
transversal power series map $ : (Jig, 9Jlg) — > (Qo; dO^o) which straightens the 
domain and puts A in the form: 

Dl + B{s, y) {fDl + D?) + T'^Ds + T'^Dy elliptic case 
A^Dl + B{s, y) {Dl - fDl) + F'^D, + F'^Dy hyperbolic case 

in the sense that the left and right sides agree to infinite order at y = 0. Here, 
B{s,y) is a transversal power series. 
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By choosing $ carefully, one can arrange that the Laplacian ($* ^A$*) i 

in the transformed coordinates has the principal term + (a22(s)^^-D| + 
en{s)yDy + boQ{s)D^. We can further put the principal part into the normal 
form 

+ |6oo (s) [y^Dl + D?] elliptic case 

(6.4) 

Dj + |6oo(s)[-D| - hyperbohc case 

6.6 Semi-classical scaling of the Laplacian 

Having straightened the domain, and hence having transferred the informa- 
tion about the boundary into the metric, we now follow the approach in the 
boundaryless case by scaling the Laplacian. Following [L, Z2], we use the 
notation N — h~^. 

We define operators Tjv, on the model space -L^(Qo) by 

. T,,f{s,y):=Nf{s,Ny) 

(6.5) 

. M^f{s,y):^e'''''fis,y) 

Definition 6.2 The semidassically scaled Laplacian is the operator on fig 
defined by 

In the straightened form, the rescaled Weyl symbol of A has the form: 

al^ ~ (a + Nf + B{s, ^y){±y^N-\a + + iVV) + K{s, ^y). (6.6) 

6.7 Conjugation to semiclassical normal form 

As in the boundaryless case, we begin by conjugating the principal term 
(which again is quadratic) 

'^I^v = + '^N^[D-s + hoo{s)i] mod N (6.7) 

into quadratic normal form by an s-dependent metaplectic conjugation on 
the transverse space to the bouncing ball orbit. Here, / denotes the quantum 
action operator: P — |(-D| + y^) in the eUiptic case and — |(-D| — y^) in 
the hyperbolic case. 
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In [Z2], we prove that there exists an 5'L ( 2, IR)- valued function aa{s) so 
that 

//(a„)*[-Dj + boois)i]ii{aa) = -D^ + fl, /^(aa)(0) = l^{ao){L) = Id, 
where a = boQ{s)ds. We then conjugate the resulting operator 

JIn := ii{aayANlJ>{aa) (6.8) 
to the semiclassical normal form 

Fjv(/)' r^N^ + N^^+ pi{i) + N-%{i) + ■■■ (6.9) 

by conjugations as in the boundaryless case, but now additionally preserving 
the boundary conditions. 

We assume for simplicity that the bouncing ball orbit is elliptic, but 
the same argument and result hold in the hyperbolic case. We again use the 
notation of (6.3): we denote by A\o the restriction of A to functions of y only. 
In the boundary case, we need to construct complex valued Weyl symbols 
Pj/2, Qj/2 SO that iterated composition with e ^ '^/"^ will successively 
remove the lower order terms in 3?jv after restriction by |o and so that the 
boundary condition is satisfied. It turns out that the boundary condition on 
the conjugating operator involves only the odd part of the real part and the 
even part of the imaginary part: 

i"j/2(o, y, n) = p-/2iL, y, n) = 0, g]/2(o, y, n) = Qj/2(^, n) = o. 

(6.10) 

We emphasize that there is no condition on the odd part Qj/2 the imaginary 
part or the even part P^^^ the real part. 

6.7.1 A Sturm-Liouville homological equation 

As before, we see what happens in the first and second steps. We first find 
Pi{s,y, Dy), Qi{s,y, Dy) so that the boundary conditions are satisfied and 
so that 

^-iV-(P+.Q), ^^^.-(P..Q). 1^ ^ ^ + ^/] + . . . }l. (6.11) 

Ij 

Expanding the exponential, we get the homological equation: 

{[Ds + yi,{P + iQ)L]+Oli}\o = Q. (6.12) 
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Taking the complete symbol of both sides we get the symbolic homological 
equation: 

i{(7+ -/,Pi +iQi} + = 0. (6.13) 

^2 2 2 

The equation may be rewritten in the form: 

d-s{P + iQ)i_{s,ra{s){y,f])) = (6.14) 
whose solution is given by 

(P + ig)i(s,r,(s)(y,7^)) = (P + iQ)i(0)-^ / -Jl.iu.y.fDdu. (6.15) 
^ ^ Jo 

We thus need to determine (P + iQ) i (0) so that the boundary conditions at 
s = and s = L are satisfied. Thus, we must solve the system 

dsPf + f {/, p?} = \o dsP? + f {/, Pf } = \o 

2 2 2 2 2 2 

(6.16) 

5.-<9t + Q 1 } = -3f^3^i lo Wi + Q 1 } = - W I, 

2 2 2 2 2 2 

with the boundary conditions 

Pf (0) = P?{L) = 0, (0) = g! (L) = 0. (6.17) 

2 2 2 2 

As in the boundaryless case, 'Ji'^{u, y,f]) = y o (ar^^{u))I is a polynomial 

2 

of degree 3 in (y, rf) in which every term is of odd degree in (^, rf) . It follows 
that Pi,Qi are also odd polynomials of degree 3. Assuming we are in the 
elliptic case, we change coordinates to the complex cotangent variables z — 
y + ifj, z — y — if] and write 



' Pf 

2 








~ 'l2m,n:m+n<3P\mni^^i^"^^^ 


+ Z'"^"), 


P? 

2 


% 


z, 








< 

2 


is 






~ Sm,n:m+n<3 ^|mn(^) 


+ Z"*^"), 


Q1 

\ 2 


is 











(6.18) 
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Then (6.16) may be rewritten in terms of these coordinates. In the eUiptic 
case, we have: 



' TsPlmni') + 'Tim 




^3^1 |o 

2 


< 




C^OJo 1 
2 




- n)q1 (s) = 


2 


j-ql (s) + if(m 


— n)qt (s) = 


2 



(6.19) 



Since m ^ n, we can ehminate p\ , q° and reduce the homological equa- 
tion to (uncoupled) second order Sturm-Liouville boundary value problems 

< -£qrnni')-ifim-^)M^ni')--^i^m^ (6-20) 

pi (0) = 0, p'i (L) = 0; ql (0) = 0, (L) = 

for the independent variable pi , gf (s). The boundary value problem 

(6.20) is always solvable unless is an eigenvalue of the operator £)| — 
[^(m — n)Y (elliptic case). The eigenfunction would have to have the form 
sin(^(m — n)s), hence a sufficient condition for solvability is that a/vr ^ Q. 
The analogous boundary problem in the hyperbolic case is always solvable. 

Thus, the conjugation eliminated the sub-principal term CRi and put our 
operator in normal form up to order N. The exponents Pi,Qi are odd 
polynomial differential operators of degree 3 with smooth coefficients defined 
in a neighborhood of [0, L] . 

6.7.2 The normal form coefficients emerge 

We carry the process forward one more step because, as in the boundaryless 

case, the even steps produce the non-trivial normal form coefficients. In out- 

1 

line, the conjugated operator in the second step, 3?^, is further conjugated 
with an exponential of the form xhe homological equation is 
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again a Sturm- Liouville boundary problem, with the new feature that diago- 
nal terms with m = n do occur in (6.19). In the elliptic diagonal term 
is a function of \z\^, hence is even under the involution z —>■ z. We write the 
diagonal terms in the form: 

{ " (6.21) 

I Qp,\z\') = Em:m<jlU''M"' 

It is impossible to solve the inhomogeneous boundary problem (6.19) for 
the diagonal terms as in the odd case with zero boundary conditions. To 
obtain a solvable system, we must add new terms /j(|2;|^) to the right side: 

iqU')--^n-j\o+m\^\') 

< (6.22) 

The right hand sides are determined by the condition that equations with 

(6.23) 

satisfies the boundary conditions. This forces 

fj{\z\') = \z?)du. (6.24) 

In the case j = 1, we have then conjugated to N^+N'^'Jl+Op^{fi{I)) + 
0{N~^). The coefficients of the polynomial /i(/) are the first Birkhoff invari- 
ants. We note that Op^{fi{I)) is a function of /, so we have conjugated to 
normal form to fourth order. We now repeat the process to complete the 
normal form. 
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6.8 Monodromy operator 

We now outline a different approach to normal forms due to lantchenko- 
Sjostrand-Zworski [SjZ, ISZ]. (See also the appendix [SjZ2]). Their goal is to 
construct the normal form of a so-called quantum monodromy operator rather 
than of the Laplacian or wave group. The normal form of the monodromy 
operator is apparently the same as the semi-classical normal form of the 
wave group that we constructed in §6.4 and §6.7. However, it is presented 
in a conceptually clearer way which works equally well in the boundary or 
boundaryless case. On the other hand, the method presented above also gives 
an algorithm for calculating the normal form, which is at present lacking in 
the approach of [ISZ]. 

The new concepts introduced in [SjZ, ISZ] following earlier one-dimensional 
definitions in works of Helffer-Sjostrand and Cohn de Verdiere- Parisse) are 
the following: 

• The monodromy operator M(A) : kerj„Q(-\/A — A) ^ keimoiVA — A) 
acting on the space kermg(-\/A. — A) of microlocal solutions of the equa- 
tion (\/A — X)u{h) = 0(A~°°) near mo, where mo is is an arbitrarily 
chosen base point of the given geodesic 7. 

• The flux norm on the microlocal solution space ker^(,(-\/A — A), with 
respect to which M(A) is unitary; 

• The Grushin reduction of the wave operator to the monodromy opera- 
tor and the expression of the wave trace as a trace of the monodromy 
operator. 

The microlocal solution space is essentially the space of quasi-modes along 
arcs of 7, so the monodromy approach is not in essence so different from the 
approach we outline above. We give a brief exposition of the key ideas, which 
are discussed in more detail in the appendix by Sjostrand-Zworski [SjZ2]. 

Assuming for simplicity that 7 projects to an embedded curve in M, we 
again work on the model space x ]R"~^ given by the normal bundle Nj 
to M. We then consider the universal cover M x M""^, fix the base point 
m = (0, 0), and consider microlocal solutions along the geodesic M x {0}. In 
the simply connected space M x M"~^, there is no obstruction to constructing 
a global quasimode along R x {0}. So we may define 

M(A) : kero(VA - A) ^ kerL(VA - A) (6.25) 
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where L — L^hy 



M(A)[iiA]o = K](L,o) 



(6.26) 



where [u\]x is the germ of the quasimode at the point x. We obtain a quasi- 
mode when (/ - M(A))[ma]o = 0. 

To make contact with the wave group, we quote the following observation 
of Sjostrand-Zworski: 

Proposition 6.3 [SjZ2j We have: 

M(A) = exp(-iL(\/A - A)) : kero(\/A - A) ^ kerL(A/A - A). 

Sketch of Proof (See [SjZ2] for more details) It is obvious that cxp(— iL(\/A— 
A)) takes microlocal solutions into themselves since it commutes with \/A. 
The rest of the proof is to verify that the complete symbol of the quasimode 
exp(— iL(-\/A — X))ux{s, y) at (0, 0) is the same as that of u\{s, y) at s — L. 
QED 

6.8.1 Flux norm 

We now describe the flux norm (•, •) on kermiV^ — A). It is motivated by 
the properties of the probability current density j of a solution ip{^^ i) of the 
Schrodinger equation in quantum mechanics (see e.g. [LL], §19). It is defined 



The integral j ■ dA over a surface measures the probability that the particle 
described by if) will cross the surface 5* in a unit of time. As is pointed out 
in [LL], liiP^ Aeii^ then j = \A\'^VS. 

The flux norm on microlocal solutions is an invariantly deflned version of 
this. We consider the microlocal solution space kcr„j(, for P = Ji^A — 1 along 
an initial arc of a closed geodesic 7. Wc let Y denote a transversal to 7 in 
M. If vy is the unit normal to Y, then 



by 



J = - 



t(^VV' -^vv-). 




More precisely and generally. 
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where we use the semi-classical notation P = —h^A — 1 with h — X~^, 
and where x is a microlocal cutoff defined in the universal cover of a small 
tube around 7, whose complete symbol equals before mg and 1 after mg 
in a somewhat smaller tube. We use the semi-classical notation to conform 
to the notation of [ISZ, SjZ] and also because their results apply to much 
more general semi-classical operators P{h). For further details we refer to 
[ISZ, SjZ]. 

Proposition 6.4 [ISZ] The monodromy operator JA, is a unitary operator 
on ker^(j(-\/A — A) with respect to the flux norm \ \ ■ \ \qf- 

6.8.2 Grushin problem 

One of the basic steps of [SjZ] is the proof of a trace formula which relates the 
trace of the wave group at 7 to the traces involving the monodromy operator. 
Formally, the trace formula says roughly that 

trp{P/h)x{P)A^y^tr I piz/h)-^logdetiI + M{z±iO))xiz)dz, (6.27) 
_i_ Jr dz 

where p e C^{Lj — e,L^ -\- e) (compare to (7.1)), and x is a cutoff to the 
sphere bundle. The precise statement is given in Theorem 1 of [SjZ]: 

trp{P/h)x{P)A^— / p{z/h)M{z,h)'-M{z,h)x{z)dz + 0{h^). 

(6.28) 

This formula is analogous to one in the case of domains with boundary, where 
the boundary replaces the transversal, see Proposition 7.1. 

To prove the formula, one sets up a microlocal Grushin problem near the 
closed trajectory 7. One begins by forming a microlocally invertible system 

(^P{z) R^{z) \ 
^ : D'{M) X 2)'(M"-^) ^ D'(M) x D'{W-^), 

R+{z) / 

where R±{z) are defined microlocally near 7. The microlocal inverse of CP(2;) 
is given by: 

( E{z) E4z) \ 
8.(z) = (6.29) 
\e_{z) E_^{z) J 
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We thus obtain the key formula 

(VA - A)-^ ~ E{X) - E+{X){I - M)~^E_(A), (6.30) 

which reduces the eigenvalue problem to the microlocal invertibility of (/ — 
M). Using this formula one proves (6.27). This equation is analogous to one 
which occurs in the Fredholm-Neumann reduction of the Dirichlet problem 
to the boundary of this type (see 7.4). 

The authors of [ISZ, SjZ, ISj] then put M(A) into quantum Birkhoff nor- 
mal form. In view of Proposition 6.3, the normal form of the monodromy 
operator is just what we called the semi-classical normal form of the wave 
group. The operation !„ thus has been interpreted as restriction to the mi- 
crolocal solution space, or as part of a Grushin reduction to the transversal. 
In [Z4, Z5] , we proved that the semi-classical normal form was a spectral in- 
variant by first turning into the homogeneous normal form, by proving that 
the homogeneous normal form was a spectral invariant and then by relating 
the semi-classical and homogeneous normal forms. The formula (6.27) elim- 
inates this latter step and directly shows that the semi-classical normal form 
is a spectral invariant by using the Grushin reduction. 

At the present time, the only application of the monodromy method to 
inverse spectral results for plane domains is for the case of bouncing ball 
orbits of analytic domains with two symmetries as in [Z2]. This docs not 
involve any new calculations since the principal symbol of the normal form 
of the monodromy operator is the classical Birkhoff normal form of CP-y and, 
as mentioned above, this normal form determines Q, when it is analytic with 
two symmetries (as proved in [GdV]). 

7 Calculation of Wave Invariants II: Balian- 
Bloch approach 

In this section, we describe a method for calculating wave invariants of 
bounded domains which does not use normal forms or parametrices, but 
rather is based on an exact formula for the Dirichlet (or Neumann) resolvent 
in terms of the 'free resolvent'. Our goal is to explain how the calculations 
of Theorem 5.2 and GoroUary 5.3 are done. 

For simplicity, we will confine ourselves to calculating wave invariants at 
a bouncing ball orbit of a bounded simply connected domain ^2 C M^. The 
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method extends with no difficulty to higher dimensions, more general metrics 
and domains (in particular, non-convex domains) and more general periodic 
orbits, but the special case is already sufficiently rich and difficult. 

The method is based on the use of a classical formula due to C. Neumann 
and I. Predholm for the resolvent Rb{z) for the boundary problem in Q 
with boundary conditions B in terms of layer potentials and the boundary 
integral operators they induce. This approach was first used by the physicists 
Balian-Bloch [BBl, BB2] in their well-known work on the Poisson relation on 
three dimensional Euclidean domains, and we will refer to it as the Balian- 
Bloch approach. The crucial advantage of this approach over the use of 
normal forms or monodromy operators is its computability. For the first 
time, one can calculate wave invariants of all orders explcitly and find out 
what information they contain. 

7.1 Wave invariants as semi-classical resolvent trace 
invariants 

First, we explain how wave trace expansions are equivalent to esolvent ex- 
pansions. We fix a non-degenerate bouncing ball orbit 7 of length L^, 
and let p G C^{L^ — e, + e) be a cutoff, equal to one on an interval 
{Ly — e/2, + e/2) which contains no other lengths in Lsp(Q) occur in its 
support. We then define the regularized resolvent by 



When 7,7~^ are the unique closed orbits of length L^, it follows from the 
Poisson relation for manifolds with boundary (see §3 and [GM, PS]) that 
the trace Trl^iRp{{k + ir)) of the regularized resolvent on L^(f2) admits a 




From the resolvent identity (e.g.) 




it follows that 




(7.2) 
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complete asymptotic expansion of the form: 



TrlnRpik + ir) ~ e^'^-^'^^~< ^{B^j + B^-ij)k-^, k ^ oo (7.3) 

j=i 

with coefficients B^ j , B^,-i j determined by the jet of at the reflection points 
of 7. The coefficients B^ j,B^-i j are thus essentially the same as the wave 
trace coefficients at the singularity t = L^. 

7.2 Reduction to the boundary 

To calculate TrlnRpik + ir) asymptotically, we use the Fredholm-Neumann 
reduction of the Dirichlet or Neumann problems in a bounded domain to the 
boundary. The key formula is similar to (6.30): 

Ra{k + ir) = Ro{k + iT)-D£{k + iT){I + N{k + iT))-^-f§£^''{k + iT), (7.4) 

where Ro{k + ir) is the free resolvent — (Aq + {k + irY)~^ on R^. Here, 
7 : H^(il) H^'^f^iydVL) is the restriction to the boundary, and D£{k + ir) 
(resp. $l{k + it)) is the double (resp. single) layer potential is the operator 
from H'{dVt) Hl^^^''^{Vt) defined by 

$£{k + iT)f{x) = GQ{k + IT, X, q)f(q)ds(q), 

(7.5) 

Diik + ZT)f{x) = J,^ £-G,{k + zr, X, q)f{q)ds{q), 

where ds{q) is the arc-length measure on where u is the interior unit 
normal to fl, and where di, = u ■ V. Also, S£*^ is its transpose (from the 
interior to the boundary), and GQ{z,x,y) is the free Green's function, i.e. 
the kernel of Ro{z). Further, 

Nik + ir)f{q) = 2 f ^Go{k + tr, q, q')f{q')ds{q') (7.6) 
Jan oi^y 

is the boundary integral operator induced by "Di. It is classical that this 
operator satisfies 

(i) N{k + iT) E'^-^dQ), 

(7.7) 

(ii) (I + N{k + ir)) : H^{dil) ^ H^{dO,) is an isomorphism. 
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The reader might now compare the formulae (7.4) and (6.30) to see how 
A^(A) is analogous to M(A). Moreover, it is used in an analogous way to 
reduce the calculation of the (distribution) trace of Rnlk + ir), formula (7.4) 
to the boundary. In doing so, we found it [ZIO] more convenient to combine 
the interior and exterior problems as follows: We write L^(M^) = L^(f2) © 
L^(Q'^) and let R^, resp. (k + ir), R^{k-\-iT) denote the Neumann resolvent 
on the exterior domain, resp. the Dirichlct resolvent on the interior domain. 
We then regard R^{k + ir) © R%{k + ir) as an operator on this space. 
The reason for combining the interior /exterior problems is that we can cycle 
around the layer potentials in (7.4) when taking the trace and simplify the 
formula to: 



Tr^2[R^\k + it) © R^{k + it) - Ro{k + ir)] = ^ logdet {^I + N{k + ir) 

(7.8)^ 

where the determinant is the usual Fredholm determinant. We refer to [ZIO] 
for a proof of this apparently well-known formula. This gives: 

Proposition 7.1 Suppose that Lj is the only length in the support of p. 
Then, 

/ p{k - A) — logdet (7 + A^(A + ir))dA ~ V B^-jk'^ , 
Je "A ^ ' 

where B^.j are the wave invariants of 7. 

7.3 Semi-classical analysis of N{k -\- ir) 

The next step is to analyse the operator N{k + ir) and the geometric series 
expansion of (/ + N{k + ir)^^. The operator N{k + ir) has the singularity 
of a homogeneous pseudodifferential operator of order —1 on the diagonal 
(in fact, it is of order —2 in dimension 2) and that is the way it is normally 
described in potential theory. However, away from the diagonal, it has a 
WKB approximation which exhibits it as a semi-classical Fourier integral 
operator with phase dgn{q, q') = \q — q'\ on dQ x dQ, the boundary distance 
function of fl. Indeed, the free Green's function in dimension two is given 
by: 

Go{k + iT,x,y) = H^'\{k + iT)\x-y\)^J^,e^(^-y'^^m'-{k + iT)')-'d^. 
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Here, '{z) is the Hankel function of index 0. It has the near diagonal and 
off-diagonal asymptotics, 



By the explicit formula wc have: 

^N{k + iT,q{cl>),q{cl>')) = d,^Go{fX, q{cl>), q{cl>')) 

= - {k + ir)H[^\k + ir\q{(l))-q{(l)')\) 
X cosZ(g((/)) - q{(f)'),iyg(^^)). 

Combining with (7.9(ii)), we sec that when 1^(0) — q{4'')\ > 1^1"*^"^ for some 
e < 1, , N{k + ir, q{(j)), q{4>')) is a semi-classical Fourier integral kernel whose 
phase is the boundary distance function. For a convex domain, the boundary 
distance function generates the biUiard map of dO, and hence we view N{k + 
ir) as a global quantization of the billiard map. For non-convex domains, 
the boundary distance function additionally generates 'ghost orbits' which in 
part exit the domain, but these only present a mild complication. We refer 
to [HZel] for discussion of these orbits and to [Z2, HZel] for further discussion 
of N{k). 

We now explain how to use Proposition (7.1) to calculate the wave trace 
coefficients B^j at a closed geodesic. When calculating these coefficients we 
first compose with a special kind of semiclassical cutoff operator x{x, k~^Dx) 
to a neighborhood of 7 on both sides of (7.4). Since we are not dealing with 
conventional Fourier integral operators, it must be proved that composition 
with such a cutoff does in fact microlocalize to 7. We will suppress this issue 
until the end. 

At least formally, we expand {I + N {k + ir))''' in a finite geometric series 
plus remainder: 

Mo 

{I+N{\+it))-^ = 5^(-l)'^ N{Xf' + N{Xf'"+\l+N{X+zT))-\ 

M=0 

(7.10) 

To calculate a given wave invariant, we need to show that, for each order k 



— 75- In \z\ as Izl — >• 0, 

ZTT II II ' 



{a) e^(^l^-2'l-'^/^)^ as |z| ^ 00. 



(7.9) 
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in the trace expansion of Corollary (7.1), there exists Mo(J) such that 



(ii) TrJ^p{k - X)N{X)^'>+\I+N{X+ir))-^N\k + ir)dX = 0{k--^-^). 

(7.11) 

We outline the method for obtaining the wave trace asymtotics at a closed 
geodesic 7 from the first term. Since N has singularities on the diagonal, one 
cannot just apply the stationary phase method to the trace. Rather one has 

to regularize the operator. We do this by separating out the tangential and 
transversal parts of N by introducing a cutoff of the form xik^"^\Q^Q'\) to the 
diagonal, where 6 > 1/2 and where x ^ C*^(M) is a cutoff to a neighborhood 
of 0. We then put 

N{k + iT) ^ No{k + iT) + Ni{k + iT), with (7.12) 

r Noik + iT, q, q') = x{k^-'\q - ?'|) N{k + ir, q, q'), 

(7.13) 

[ Ni{k + ir, q, q') = (1 - x(A;^"^|g - g'|)) N{k + ir, g, q'). 

The term A^i is a semiclassical Fourier integral kernel quantizing the billiard 
map, while A'^o behaves like an Airy operator close to the diagonal with the 
singularity of a homogeneous pseudodifferential operator on the diagonal. 

Now consider the powers N{k + it)^ in the first term (i) of (7.11). We 
write 

(Aro + Ar,)^= N,ii)oN,(2)0---oN,^M). (7.14) 

f7:{l,..,M}^{0,l} 

We regularize A^*^ by eliminating the factors of A"o from each of these terms. 
This is obviously not possible for the term A"^^ but it is possible for the 
other terms. By explicitly writing out the composition in terms of Hankel 
functions and using the basic identities for these special functions, we prove 
that A^o°-^i°Xo(^ + ^''', 01, ^2)) is a semiclassical Fourier integral operator on 
dQ of order —1 associated to the billiard map. Thus, composition with A^o 
lowers the order. The remaining terms N^^ , when composed with a cut-off 
to 7, do not contribute asymptotically to the trace. 
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The successive removal of the factors of A^"o thus gives a semi-classical 
quantization of the billiard map near 7. We then calculate the traces of 
each term by the stationary phase method and obtain the result stated in 
(5.2). The terms displayed there with the maximum number of derivatives of 
the defining function of dQ come only from the Nf^ terms. Thus, although 
the removal of the A^o factors does change the amphtude of the Ni factors 
and does contribute to the wave trace, it turns out to be negligeable in the 
inverse spectral problem. This is one of the principal virtues of the Balian- 
Bloch approach. 

We will explain how the explicit formulae for the wave invariants are 
derived in the next section. 

7.3.1 Remainder estimate 

We now address some sketchy remarks to the remainder estimate, in the 
hopes of convincing the reader that it is plausible to expect the remainder 
trace to be small. The ability to insert a microlocal cutoff to 7 is crucial here. 
The main obstacle to the remainder is that the norm of N{k + ir) fails to 
decrease with increasing r due to the Airy part associated to creeping rays. 
This again is a difference to the monodromy operator. 

To obtain a small remainder we set the spectral parameter in equal to 
k + ir log k. The presence of the log k in the imaginary part changes the wave 
trace expansions by k"'--^^^^ , but this does not hurt the expansions such the 
remainder estimates will be of lower order. We then estimate the remainder 

Tr I p{k - X)N{X)^°+\l+N{X+iT))-^ Xj N'{k + iT)dX 
Jr 

by applying the Schwarz inequality for the Hilbert-Schmidt inner product, 
and using the relation (/+iV(A+ir))"^ to the Poisson kernel to estimate this 
factor. These estimates leave a trace of N^N*^ for fixed M, microlocahzed 
to 7. We regularize these traces as above and obtain osillatory integrals 
whose phases have critical points corresponding to M-link closed circuits 
which being at some point q, bounce along the boundary until q' and then 
return to q by traversing the links in reverse order. However, the cutoffs to 
7 force the links in critical paths to point in the direction of 7 and hence to 
be of length roughly ML^. The imaginary part it log k of the semiclassical 
parameter then contributes a damping factor of ^-'^^^-r'^^sk each link. 
The links correspond to the A^i factors. Thus, for each string, we have one 
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k ^ for each Nq factor and one e '^'^^t^°sk fQj, gg^^j^ factor. For sufficiently 
large r these combine to give a factor of for any prescribed R. 

7.3.2 N{X) versus M(A) 

We digress momentarily to compare A^(A) and the monodromy operator 
M(A). Both operators arise as a reduction of the wave group to the boundary 
(or a transversal) and are quantizations of the billiard map, N{X) globally 
and M(A) microlocally at a closed orbit 7. (An earher reduction also occurs 
in [MM, P, P3] but in a rather different way. ) In applications to boundary 
problems, they are quite similar, as evidence by the comparisons of (6.27) 
with Proposition 7.1 and of (6.30) with (7.4). 

But there do exist rather important differences. The operator iV(A) is 
not a standard Fourier integral operator on dQ, while M(A) is one. The 
diagonal singularities of A^(A) require a complicated rcgularization procedure. 
Moreover it is a global invariant of dQ, not a microlocal one at 7; On the 
other hand, A^(A) is just the restriction of a canonically defined free Green's 
kernel to dfl x dQ, and thus is an elementary and computable object. By 
comparison, M(A) must be constructed by some kind of parametrix method. 
It was precisely the complexity of micrlocal parametrices for the wave group, 
even at periodic reflecting rays, which motivated our turning to the Balian- 
Bloch approach. 

7.4 Stationary phase expansion 

After regularizing the traces, we end up with oscillatory integrals in the 
standard sense and obtain expansions by applying stationary phase. The 
amplitudes and phases are canonical, since we began with canonical ampli- 
tudes and phases and since the rcgularization procedure is essentially the 
same for all domains. So most of the complexity of the expansion is due to 
the stationary phase method. 

The key point in inverse spectral theory is to identify data in the station- 
ary phase term of order k^^ which represents 'new data' not contained in 
the previous terms. Terms of the coefficient of k~^ which contain the maxi- 
mum number of derivatives of the phase are the most important ones. Thus 
we face the combinatorial problem of locating such terms in the stationary 
phase expansion. The Feynman diagram method of assigning labelled graphs 
to each term in the expansion proves to be very effective for this purpose. 
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Consider a general oscillatory integral Zj. = a{x)e^'^^^^^dx where a e 
C^(]R") and where S has a unique critical point in suppa at 0. Let us write 
H for the Hessian of S at 0. The stationary phase expansion takes the form: 



Here, the sum runs over the set Sv,i of labelled graphs (F, £) with V 
closed vertices of valency > 3 (each corresponding to the phase), with one 
open vertex (corresponding to the amplitude), and with / edges. Further, 
the graph F' is defined to be F minus the open vertex, and xr = V — I equals 
its Euler characteristic. We note that there are only finitely many graphs for 
each X because the valency condition forces I > 3/2V. Thus, V < 2j, I < 3j. 

The function £ 'labels' each end of each edge of F with an index j G 
{1, . . . , n}. Also, S{r) denotes the order of the automorphism group of F, and 
/^(F) denotes the 'Feynman amplitude' associated to (F,£). By definition, 
Ie{T) is obtained by the following rule: To each edge with end labels j, k one 
assigns a factor of ^h^'^ where — [y^). To each closed vertex one assigns 

a factor of ifc ^^i^ ^^^^ where v is the valency of the vertex and ii. . . ,1^, at 
the index labels of the edge ends incident on the vertex. To the open vertex, 
one assigns the factor g^fi " g^i,. ; where z/ is its valence. Then /^(F) is the 
product of all these factors. To the empty graph one assigns the amplitude 
1. In summing over (F, t) with a fixed graph F, one sums the product of all 
the factors as the indices run over {1, . . . , n}. 



An analysis of the diagrams and amplitudes shows that the jth even Taylor 
coefficients f^^\Qi) of the boundary defining functions appear first in the 
k^^^^ term. 

When the domain has one symmetry axis, which we visualize as an 
up/down symmetry, the terms with this data have the form 




7.4.1 The data /±^(0) 



2rL(/ii^)V^''nO) + • • • , 



where • • • refers to terms with <2j — 1 derivatives. 
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When the domain has two symmetries, a left / right symmetry in addition 
to an up/down symmetry, the odd Taylor coefficients vanish and we see 
immediately that the even Taylor coefficients can be determined inductively 
from the wave trace invariants. This gives a new proof that analytic domains 
with two symmetries can be determined from the wave trace invariants at a 
bouncing ball orbit which is one of the symmetry axes. 

This does not quite prove that such domains are spectrally determined 
among other analytic domains with two symmetries, since the length of the 
bouncing ball orbit must be known in order to obtain the wave invariants. 
This length is a spectral invariant if the domains are additionally convex 
[Gh]. For non-convex domains we needed to add as an assumption that the 
bouncing ball symmetry axis had a fixed length L in [Zl] . 

7.4.2 The data /4^^~^\o) 

When we do not assume a left /right symmetry, the odd Taylor coefficients 
are non-zero in general, and the problem arises whether there is sufficient 
information in the wave invariants to determine all of the even and odd 
Taylor coefficients of the boundary defining function (or curvature function) 
of a domain with one symmetry. We assume the axis of symmetry is a 
bouncing ball orbit whose orientation is reversed by the symmetry. 

An analysis of the diagrams and amplitudes for the odd Taylor coeffi- 
cients ^■'(0) show that they appear first in the term of order fc^-'^^. It 
turns out that five diagrams contain this data, but the amplitudes of three 
automatically vanish. They two amplitudes have the following forms: 

• (i) {hfy-^h'!^hF^f^^-^\Q)f^\Q). 

• (ii) {hfy-'m'f^^'-'\Q)f^'\Q). 

To decouple them, we need to analyze the behavior of power sums of 
columns the Hessian matrix elements. In [Z5] we proved that cubic column 
sums are linearly independent from linear ones as r ^ oo. 

7.5 Positive results for analytic domains and metrics 

We now review the proof in [Z5] that the Taylor coefficients /^■'"^(O), /^■'(O) 
can be determined from the wave invariants B^,] as r varies over r = 
1,2,3,.... 
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It suffices to separately determine the two terms 



A . ^ (7.15) 

and [ElUhll)y^'\0)f('^-'\Q). 

It is easy to see that the terms decouple as r varies if and only if the cu- 
bic sums X]gLi(^2r)^ are non-constant in r = 1,2,3, By the explicit 

calculation in [Z2], we have: 

2^2r 1 
^ki,k2=0 (cosha/2+cos ^)(cosha/2K;os ^ ) (cos a/2tcosh '^''-^-^^i)^ )' 

It is obvious that the sum is strictly increasing as r varies over even integers. 

This independence of the linear and cubic sums permits us to use an 
inductive argument when the domain has the given symmetry. From the 
J = term we determine f"{0). Indeed, (1 — L/*^^^(0) = cos{h)a/2 and a is 
a wave trace invariant. From the j = 2 term we recover /'^(O), /^(O). The 
induction hypothesis is then that the Taylor polynomial of / of degree 2j — 2 
has been recovered by the j — 1st stage. By the decoupling argument we can 
determine /^-'(O), P^~^{0). (Strictly speaking, there is the minor annoyance 
of the factor of /*-^^(0), which is resolved in [Z5]). 



7.5.1 Final comments on the Balian-BIoch invariants 

It would of course be desirable to remove all the symmetry assumptions, if 
possible. Thus one would need to recover the Taylor coefficients /i^^^O), f±^~^\o) 
when /_)_ 7^ /_ from the wave invariants. The problem is that the induction 
does not work because the term i?2r(j^-'^^/(0)), which we do not need to 
know in the symmetric case, might have a different dependence on the Tay- 
lor coefficients of f± than the 'principal term', i.e. the one with the highest 
derivatives. The latter is basically a sum of the coefficients of /± while 
R2r{i^^~^ f{^)) could involve any symmetric polynomial in these coefficients. 

(2) ('2) 

At the present time, we do not even know whether the set {/| ^0), /- '(0)} 
of second derivatives can be determined from the second wave coefficient. 

Some hope has been provided by computer calculations of C. Hillar [Hil]. 
The results suggest that Hessian column power sums with different powers 
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are linearly independent as r — > oo. This would give quite a large supply of 
Taylor coefficients. 

8 Surfaces and domains with integrable dy- 
namics 

As mentioned above, Birkhoff normal forms at a periodic orbit (on both the 
classical and quantum level) are approximations to the Hamiltonian by an 
integrable one. When the Hamiltonian is completely integrable, its dynamics 
and Birkhoff normal forms are very special and the metrics might be rather 
spectrally rigid, at least in low dimensions. This is the case with flat metrics 
[Ku3]. In this section, we consider some model inverse spectral problems on 
two-dimensional surfaces or domains with integrable dynamics. 

Let us recall that the geodesic flow of an n-dimensional Riemmanian 
manifold (M, g) is completely integrable if it commutes with a Hamiltonian 
action of on T*M - (n = dimM). That is, the metric Hamiltonian 
satisfles 

{l^lg^Pj} = = {PhPj}, i,j = l,...,n 

where pj : T*M — —>■ M. are homogeneous of degree one and are independent 
in the sense that 

dpi A dp2 A • • • A dpn 7^ on a dense open subset U cT*M. 

The orbits of an R" action give a (usually singular) fohation of S*M (the 
unit sphere bundle for the metric), called the Liouville foliation, by afflne 
manifolds of the form M."- ■ (x,^) = R"^/r where F is the isotropy group at 
(x, ^). If it is a lattice of full rank, the orbit is a torus T" of dimension n. If 
it has less than full rank, the orbit type is T'^ x R*^"*^. The isotropy group 
might have the form R'^ x Z"~'^ in which case the orbit becomes a singular 
torus of dimension n — k. The orbits which contain periodic geodesies are 
sometimes called 'periodic tori', and they furnish the components 7 in the 
trace formula. 

In the case of bounded plane domains, complete integrability of the bil- 
liard map P means that B*dQ is foliated by invariant curves. The well-known 
conjecture of Birkhoff is that ellipses are the only example of compact smooth 
Euclidean plane domains with integrable billiard map. There are more ex- 
amples if curved metrics are allowed, see Popov- Topalov [PT] . 
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8.1 Trace formulae for integrable systems 

So far, we have mainly considered the trace of the wave group around non- 
degenerate periodic orbits. For integrable systems, the periodic orbits usually 
come in families filling out invariant tori. We now consider the appropriate 

notions of non-dcgcneracy in this context. 

We will always assume that the closed geodesies come in clean families in 
the following sense: 

Definition 8.1 ^ metric g on a compact manifold M will be said to have 
a simple clean length spectrum if the length function Lg on the loop space 
Map{S^, M) is a Bott-Morse function which takes distinct values at distinct 
components of its critical set, Crit{Lg). 

The term Bott-Morse means that each component of Crit{Lg) is a mani- 
fold, whose tangent space is the kernel of dLg. Equivalently, each component 
is a clean fixed point set for . One needs the clean (Bott-Morse) condition 
to get a nice wave trace expansion and one needs the simple length spectrum 
condition to determine geometric information from the expansion. 

Under this assumption, the trace of its wave group has the form: 

Tre'*VA;^eo(i) + J^e^W (8.1) 

where the singular term eo(t) = C„ Vol{M, (7) (t + iO)""" -|- . . . at t = is the 
same as in the non-degenerate case, where {T} runs over the critical point 
components of L^, and where 

e-j = c-j^dT {t - Lt+iO) -"^/^ + cr,d^-i{t - W+iO) ^-<'^'^+^) + . . . , d<j^d\m7. 

Here, is the dimension of the symplectic cone formed by the family of 
closed geodesies within T*M and L'j is the common length of the closed 
geodesic in T. For instance, in the non-degenerate case, T = K+7 is the 
symplectic cone generated by 7 C 5"*^^ and d<j — 2. 

8.2 Spectral determination of simple surfaces of revo- 
lution 

The simplest example of what can be done with integrable systems is given 
by 'simple' analytic surfaces of revolution. We first review the proof in [Z2] 
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that they are spectrally determined among other simple analytic surfaces of 
revolution and then sketch a result observed independently by the author and 
G. Forni and by K.F. Siburg regarding their spectral determination among 
all metrics on S"^. 

The precise class of metrics we consider are those metrics g on S'^ which 
belong to the class Jl* of real analytic, rotationally invariant metrics on 
with simple length spectrum in the above sense and satisfying the following 
'simplicity' condition 

• g ^ dr"^ + a{rfde'^ 

• 3!ro : a'(ro) = 0; 

• The Poincare map CPq ol r — Tq is elliptic of twist type 

Convex analytic surfaces of revolution are examples, but of course there 
are others. The unique isolated closed geodesic (at distance Tq) is an eUiptic 
orbit. Peanut surfaces are obviously non-simple, since they possess three 
isolated closed geodesies of which one is hyperbolic. The standard round 
metric and other Zoll metrics of revolution are non-simple since the Poincare 
maps are not twist maps. 

The following results shows that one can solve the inverse spectral prob- 
lem in the class of analytic simple surfaces of revolution. 

Theorem 8.2 ([Z2]) Suppose that gi, g2 are two real analytic metrics on S'^ 
such that {S^,gi) are simple surfaces of revolution with simple length spectra. 
Then Sp(Ag^) — Sp(Ag^) implies gi — g2- 

The proof is based on quantum Birkhoff normal forms for the Laplacian 
A. But the special feature of simple surfaces of revolution is that there exists 
a global Birkhoff normal form as well as local ones around the critical closed 
orbit or the invariant tori. This is because A is a toric integrable Laplacian 
in the following sense: there exist commuting first order pseudo-differential 
operators /i,i2 such that: 

• the joint spectrum is integral, i.e. Sp{3) C n F -|- {/j,} where F is the 
cone I2 > \h\ in M^. 

• The square root of A is a first order polyhomogenous function y/A — 
H{Ii,l2) of the action operators. 



79 



By polyhomogeneous, we mean that H has an asymptotic expansion in 
homogeneous functions of the form: 

Hr^Hi + Ho + H_^ + ..., Hj{rl) = r^Hj{I). 

The principal symbols Ij of the Jj's generate a classical Hamiltonian torus 
action on T*S'^ — 0. Analysis of the normal form shows that Hg = 0. 
It follows that 

Spiy/A'g) = {H{N + //) : TV e n Tj, 

where the eigenvalues have expansions 

Aiv ~ Hi{N + + H_i{N + //) + .... 

Theorem 8.3 ([Z2]) Let {S'^,g) be an analytic simple surface of revolution 
with simple length spectrum. Then the normal form if (^1,^2) is a spectral 
invariant. 

The normal form and the proof are very different from the non-degenerate 
case in [G, Z3, Z4], although the philosophy of the proof is similar. 

To complete the proof of Theorem 8.2, we need to show that H determines 
a metric in IR. As in the bounded domain case outline above, the crucial point 
is to calculate the normal form invariants. It turns out to be sufficient to 
calculate Hi = H and in terms of the metric (i.e. in terms of a{r)) and 
then to invert the expressions to determine a(r). 

The method given in [Z2] for calculating H and was to study the 
spectral asymptotics of -\/A = H{Ii, I2) along 'rays of representations' of the 
quantum torus action, i.e. along multiples of a given lattice point {no,ko). 
The lattice points {uo, kg + |) parametrize tori Tn^^^o satisfying so-called 
Bohr-Sommerfeld quantization conditions, which imply that one can con- 
struct associated joint eigenfunctions (f)n„,ko of (-^i; -^2) by the WKB method. 
This is reminiscent of the quasi-mode method mentioned above around non- 
degenerate closed geodesies, but there we studied the conjugation to normal 
form rather than the asymptotics of quasi-modes. Here, the existence of a 
global torus action makes the quasi-modes easier to study, and indeed the 
{<Pn,k} are actually modes (eigenfunctions) of A with complete asymptotic 
expansions along rays. By studying the eigenvalue problem as | (ri, A;) | — > 00 
we determine the H_/s. 
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Once it is known that the global BirkhofF normal for if is a spectral 
invariant, it follows that for each Uq G Z, the function H{no, I2) is a known 
function of the variable I := l2- Its principal symbol Hn^[I) := H\{no-,I) is 
then a known function and its inverse function 



l-r+{E) I 7 

.XE)^ Je- — 

Jr-{E) V 



-dr 

is also known. Here, r±{E) are the upper and lower values of the radius of 
the projection of the torus of energy E and angular momentum no to S"^. We 
may write the integral in the form 



where /i is the distribution function fi{x) :— \{r : < of ^, with 
the Lebesgue measure. This Abel transform is invertible and hence 



r: ^ 



and therefore 

J{x):= J2 



\a'{r)\ 



r:a(^r)=x 

are spectral invariants. By the simplicity assumption on a, there are just two 
solutions of a(r) = x] the smaller will be written r_(x) and the larger, r+(x). 
Thus, the function 

/ , , , 1 1 

\a'{r-{x))\ \a'{r^{x))\ 

is a spectral invariant. 

By studying we find in a somewhat similar way that 

K{x) = \a'{r_{x))\ + \a'{r+{x))\ 

is a spectral invariant. It follows that we can determine a'(r+(a;)) and 
a'(r_(x)). 

Since both metrics gi and g2 are assumed to belong to 3^*, they are 
determined by their respective functions %(r). We conclude that ai — 02 
and hence that Qi— 92- 
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8.3 Unconditional spectral determination 

We now ask whether we can remove the assumption that g2 G 31* in Theorem 
8.2? The question is, if SpecAg = SpecA/j and g e Oi*, then is h e An 
affirmative answer would give a large class of metrics which are spectrally 
determined. To the author's knowledge, the only metric on S'^ known to be 
spectrally determined is the canonical round one. 

We cannot answer this question, but we will give a partial result sug- 
gesting that it is true. The following theorem was stated in [Z7] and was 
worked out in a conversation with G. Forni in 1997. A similar but some- 
what stronger conclusion was drawn by K. F. Siburg (Theorem 5.2 of [S2]), 
under the stronger hypothesis of an isospectral (or length-spectral) defor- 
mation. After sketching our proof, we will also sketch his, which has other 
applications. 

Theorem 8.4 Let g E 01* and suppose that h is any metric on with 
simple clean length spectrum for which SpecAg = SpecA^. Then h has the 
following properties: 

• (i) It has just one isolated non- degenerate closed geodesic (up to 
orientation); all other closed geodesies come in one-parameter families 
lying on invariant tori in S^^S'^; 

• (a) The Birkhoff normal form of G\ at is identical to that of G*g at 
its unique non- degenerate closed orbit. Hence it is convergent. 

• (Hi) the geodesic flow G\ of h is C^- integrable. That is, S^S"^ has a 

-foliation by 2-tori invariant under G\. 

If we knew in (ii) that the Birkhoff transformation conjugating to 
its Birkhoff normal form was convergent, then it would follow that is 
completely integrable with global action-angle variables, and that it would 
commute with a Hamiltonian torus action. We conjecture that this is the 
case. Statement (iii) shows that it is at least integrable in the sense. At the 
present time, metrics on S'^ whose geodesic flows commute with Hamiltonian 
torus actions have not been classified. For a result which classifies such 
metrics on the torus (they must be fiat), see [LS]. 
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8.3.1 Wave invariants for g e Ji* 

The proof is based on a study of the wave trace formula in this setting. 

Proposition 8.5 Suppose that g & 01* . Then the trace of its wave group 
has the form: 

Tre^t^^ = eo{t) + e,^{t) + J] ej{t) 

T 

where 

eo{t) = Cnarea{M, g){t + iO)~" + . . . 
is singular only at t — 0, where 

e^g{t) = Cj{t — + iO)~^ + a^o log(t — + iO) + . . . 

and where 

eT=CT(t-LT+iO)-^/2 + .... 

Suppose now that h is any other metric with Spec — SpecA^ and 
with simple length spectrum. Then the wave trace of h has precisely the 
same singularities as the wave trace of g. Since there is only one singu- 
larity of the order {t + iO)~^, there can exist only one non-degenerate closed 
geodesic, proving (i). By Guillemin's inverse result, the BirkhofF normal form 
of the metric and Poincare map for 7/1 is the same as for 7^. In particular, 
the Poincare map CP^ of 7/1 is elliptic of twist type. From the fact that all 
other critical components have the singularity of a three-dimensional cone, it 
follows that in 5*^5*^ the other closed geodesies come in one-parameter fam- 
ilies. They weep out a surface foliated by circles, which can only be a two 
dimensional torus, proving (ii) 

It is the third statement (iii) which requires a new idea. So far we only 
know that periodic orbits lie on invariant tori, but we do not know what lies 
between these tori. Aubry-Mather theory will now close the gaps. 

We recall that Aubry-Mather theory is concerned with an area-preserving 
diffeomorphism of an annulus A ^ x {a,b) [KH, MF, Ka, M, SI, S2]. 
Let (f) denote a lift to R x (a, b) with + l,y) — (j){x, y) + (1, 0). The map 
(j) is called a monotone twist mapping if it preserves the orientation of A, 
if it preserves the boundary components and if the lift 4>{xQ,yQ) = {xi,yi) 
satisfies: 

• The twist condition: |^ > 0; 
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• The exactness condition: yidxi — i/odxo — dh{xo,Xi). 

If a, b are finite, (p extends continuously to the boundary as a pair of 
'rotations': 

(f){x, a) = {x + a;_, a), ; 4>{x, b) = {x + a;+, b). 

The interval (a;_,a;+) is called the twist interval of (p. 

Let {{xijUi}} be an orbit of (f). Its rotation number is defined to be 

|i|-»oo i 

A curve C G A is called an invariant circle if it is an invariant set which 
is homeomorphic to the circle and which separates boundary components. 
According to Birkhoff 's invariant circle theorem, an invariant circle is a Lip- 
schitz graph over the factor of A. Any invariant circle has a well-defined 
rotation number (the common rotation number of orbits in the circle) and 
the rotation number belongs to the twist interval. 

An orbit {(xj, yi)} is determined by the sequence {xi} of its x-coordinates. 
It is called minimal if every finite segment is action-minimizing with fixed 
endpoints: 

n— 1 n— 1 

^ h{Xi, Xi+i) < ^ h{^i, ^j+i), V(^fc, . . . , Cn) with = Xk, = Xn- 

i=k i=k 

The corresponding orbit orbit {xi,yi) is called a minimal orbit. The Aubrey- 
Mather theorem states (cf. [Ka], Theorem 1;) : 

A monotone twist map possesses minimal orbits for each rotation number 
uo G (c<Ji,cj_|_) in its twist interval. Every minimal orbit lies on a Lipschitz 
graph over the x-axis. For each rational rotation number uj = ^ , there 
exists a periodic minimal orbit of rotation number | . When cu is irrational, 
there exists either an invariant circle with rotation number cu, or an invariant 
Cantor set E. 

The theorem also describes three possible orbit types in both the rational 

or irrational case. 

We now return to our problem on simple surfaces of revolution. We fix 
local transverse discs (Poincare sections) 5*^, resp. Sh to the geodesic flows 
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G*, resp. G\ at the orbits 7^, resp. 7^. Concretely, the transversals can 
be taken to be small variations of 7^, resp. 7^^ moved up and down a small 
orthogonal geodesic arc to 7^, resp. 7^. Since 7^ is non-degenerate elliptic, its 
Poincare map defines area-preserving map of the symplectic disc Sg with 
a non-degenerate elliptic fixed point corresponding to 7^. Since the Birkhoff 
normal forms of and are the same, 7^ defines the same kind of map 
of S]i. To obtain a twist map of an annulus, we puncture out the fixed point 
of resp. J/j. We define the rotation angle ujq corresponding to this orbit 
by continuity from nearby orbits, which can be read off from the Birkhoff 
normal form. Indeed, the Poincare maps have the form 



X 



y 



cos 27:6 — sin 27r6 \ 






sin 27:9 cos 27:9 J 


M 


V J 



+ 0{\x\' + \y\'), 



with 6 ujo + (3{x'^ + y"^) + ■ ■ ■ , as a;^ -|- — > 0. For further discussion of 
area- preserving maps around elliptic fixed points and twist maps, see [S2]. 

The foliation of S*Sg\^h by 2-tori intersects Sg in a foliation by invariant 
circles converging to the fixed point. Circles with rational rotation numbers 
contain only periodic orbits, and conversely all periodic orbits belong to 
invariant circles with rational rotation numbers. 

In the case of /i, we know that once 7^ is punctured out, all periodic orbits 
come in two-tori which project to invariant curves in Sh which are diffeomor- 
phic to circles. The Birkhoff normal form shows that they are invariant circles 
in the sense that they are also homotopically non-trivial in the punctured Sh 
(indeed, the action-angle variables are essentially polar coordinates on this 
disc). 

Let C <Z Shhe an invariant circle for J/j. Set: 

• I{C) — area enclosed by C. 

• Ci — invariant circle enclosing an area /. (It is clearly unique). 

• LUi = rotation number of T'hIci- 

Since CP/j is a twist map, the rotation number cui is a monotone increasing 
function of /. 

Let /+ denote the area of 5*^ with respect to the symplectic form. By 
shrinking S^ we may assume that the boundary of Sh is a periodic circle Ci_^_. 
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Since the origin is a fixed point of CP^ and since tlie rotation number 
is increasing, it is clear tliat the set of rotation numbers lie in the interval 
[0,0'/^]. It follows by the Aubry-Mather theorem that every rational number 
p/q e [0,a;/_^] is the rotation number of a periodic circle Ci C Sh- 

Lemma 8.6 Let a e [0, Then there exists an invariant circle for of 
rotation number a. 

Proof This follows from Corollary 6.1 of [Ka], which shows that as soon 
as Birkhoff periodic points of all rational rotation numbers are constructed, 
then there exist orbits of each irrational rotation number which are dense in 
an invariant circle. Alternatively one could let Pn/qn — ^ « and let Cp^/g^ be 
the corresponding periodic circles. Each Cp^/g^ is a Lipschitz circle and the 
sequence of these circles tends monotonically to a limit circle. It is Lipschitz 
and its rotation number is a. 

To complete the proof, we make the observation: 

Lemma 8.7 Sh is foliated in the C° sense by invariant circles for CP/j. 

Proof: If not there exists an annulus A G Sh with boundary consisting 
of two invariant circles and containing no invariant circles in its interior. 
But by the Aubrey- Mather theorem, there must exist a periodic point in A. 
Since the periodic points come in circles, there must exist a periodic circle, 
contradicting the non-existence of invariant circles in A. QED 
This completes the proof of Theorem 8.4. 

8.3.2 Isospectral class of an ellipse 

Instead of surfaces in D^*, one can apply this reasoning to the Dirichlet (or 

2 2 

Neumann) problem for an ellipse Ea,b = {{x,y) '■ fs- + = !}• It is well- 
known that eUipses have integrable billiards. 

The eUipse has three distinguished periodic billiard orbits: 

• The bouncing ball orbit along the minor axis, which is a non-degenerate 
elliptic orbit; 

• The bouncing ball orbit along the major axis, which is a non-degenerate 
hyperbolic orbit; 

• Its boundary. 
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All other periodic orbits come in one-parameter families. The existence of 
a hyperbolic orbit means that the billiard flow is very different from geodesic 
flows of metrics in "Jl*. 

The wave trace formula shows that any domain fl with Spec{fl) = Spec{Ea,h) 
has precisely one isolated elliptic orbit, one isolated hyperbolic orbit. The 
accumulation points in the length spectrum must be multiples of the perime- 
ter of the domain (a spectral invariant), so the boundary must be a closed 
geodesic as well. 

One can apply the twist map theory either to the boundary orbit or to 
the unique non-degenerate elliptic orbit with isolated length in the length 
spectrum. The argument above shows that there exists a C° foliation by 
invariant circles at least near these two orbits. 

8.4 Marked length spectral rigidity of domains 

We now review a result due to K.F. Siburg [SI, S2] on isospectral deforma- 
tions of integrable systems which has interesting applications to metrics in 
%* and also to bounded plane domains. 

The key invariant is the mean minimal action 



of a twist map 0, which associates to a rotation number u in the 'twist 
interval' the mean action 



of a minimal orbit {qi,r)i) of (f) of rotation number cu. It is a strictly convex 
function which is differentiable at all irrational numbers. If uj = p/q. then a 
is differentiable at u if and only if there exists an invariant circle of rotation 
number p/q consisting entirely of periodic minimal orbits. If a monotone 
twist map possesses an invariant circle of rotation number a;, then every orbit 
on the circle is minimal ([MF], Theorem 17.4). In the case of a bounded plane 
domain, h{q, q') = —\q — q'\. 

It is observed by K. F. Siburg [SI] (Theorem 4.1) that the marked length 
spectrum is essentially the same invariant as the mean minimial action. The 
mean minimal action is therefore an isospectral deformation invariant. He 



(8.2) 




(8.3) 



i=-N 
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used this to give a proof that there cannot exist isospectral deformations 
within 01* (or within more general classes of deformations, see [S2], Theorem 
4.5). Indeed, the mean minimal actions a would all be the same. Hence if that 
cto of the original surface is differentiable at all rationals, so are they all. But 
this implies that they all have invariant circles of rational rotation numbers. 
By taking limits, one obtains invariant circles of all rotation numbers. (We 
remark that the trace formula already shows that there existed invariant 
circles of all rational rotation numbers). 

Siburg further connects the mean minimal to the Melrose-Marvizi invari- 
ants: Let a* : [—1,1] — > M denote its convex conjugate of a. Then Siburg 
shows (loc. cit. p. 300) that the Melrose-Marvizi invariants are algebraically 
equivalent to the Taylor coefficients of (a*)^/^ at —1. 

It appears that the only explicitly known mean minimal action is that of 
the disc D, where it is given by a{u!) = ^sinTra; ([SI]); a is only smooth 
when fl — D ([SI], Theorem 4.6) It is hkely that it is computable in the 
case of an ellipse (perhaps in terms of elliptic functions). Perhaps it can be 
proved that ellipses are the unique domains with these particular minimal 
action functions (which could be simpler than the long outstanding problem 
of proving that they are the unique integrable billiard systems). If so, this 
would prove that ellipses are spectrally rigid. 

Problem 8.1 Is the map from curvature functions n of convex plane do- 
mains to the mean minimal action a of the associated convex domain injec- 
tive or finitely many to one? at least near ellipses or under some additional 
analyticity or discrete symmetry condition? 
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